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1.1  OBJECTIVES 

~ — The  objective  of  the  Elast  OverPressure  (BOP)  program  is  to  understand 
the  mechanisms  by  which  biological  organa  are  injured  so  that  a  Damage  Risk 
Criteria  (DRC)  for  humans  can  be  defined  in  terms  of  blast  parameters.  In 
addition  to  the  various  experimental  projects,  mathematical  models  have  been 
developed  to  understand  the  underlying  mechanisms,  correlate  mechanical 
response  and  observed  injury,  and  to  provide  a  way  to  extrapolate  animal 
experimental  studies  tc  the  prediction  of  human  exposure  safety  guidelines. 

This  report  documents  the  construction  of  mathematical  models  to  study 
animal  response  to  external  blast.  Based  on  sheep  anatomy  a  simplified  16- 
element  thiee-dimenslcnal  model  was  constructed  to  capture  the  gross  body 
respouse.  A  refined  two-dimensional  model  was  then  constructed  to  follow  wave 
propagation  inside  the  thorax  and  to  identify  locations  subjected  to  the 
greatest  stress.  <L  ^ 

1.2  1ACKGRODHD 

Extensive  experiments  of  animal  exposure  to  eir  and  underwater  blast 
wavee  have  been  carried  out  by  the  Lovelace  Foundation  of  Albuquerque,  New 
Mexico  in  the  past  30  years  (10,  15,  16,  17].  The  study  of  blast  effect  on 
biological  bodies  has  also  been  advanced  by  the  Swedish  group  led  by  Clemedson 
and  Jbnsson  [3,  4,  E,  18,  19]. 

In  support  of  these  experiments,  Bowen  et  a'i.  [15]  in  1965  proposed  a 
mechanical  model  of  the  lung  for  studies  in  blast  biology.  In  this  model  the 
lung  is  considered  to  be  an  air-filled  cavity  with  the  chest  and  abdominal 
walls  acting  as  pistons  which  move  under  the  load  of  the  blast  wave  and 
compress  the  air.  Physically,  the  inertia,  elasticity  and  damping  effects  from 
various  organs  and  parts  of  the  body  are  lumped  into  mass,  spring,  and  dashpot 
representations.  Calculated  pressure  in  the  air  cavity  represents  the  average 
pressure  in  the  whole  lung  at  any  given  time  when  the  body  is  under  blast 


exposure.  This  idealised  model,  however,  ignores  the  wave  propagation  charac¬ 
teristics  Inside  the  body  thorax  and  its  possible  links  to  the  resulting  body 
injuries. 

Clinical  and  experimental  examination  of  blast  Injured  lung,  however, 
suggest  wave  dynamics  may  play  an  important  role  in  causing  damage.  With  the 
differences  in  wave  speeds  between  the  lung  and  the  other  organs  [6].  wave 
propagation  is  considered  an  important  phenomenon  in  the  body  during  blast 
exposure.  Reflection  and  transmission  of  incident  wave  can  be  very  different 
due  to  the  mismatched  Impedance  of  the  propagation  media.  For  a  given 
compression  loading,  spatial  and  temporal  variations  of  the  overpressure 
response  at  different  regions  of  the  lung  are  expected. 

With  this  background,  a  finite  element  approach  was  taken  in  the  work 
described  in  this  report.  With  it  we  are  able  to  see  various  wave  propagation 
characteristics.  e.g.f  focusing,  reflection,  inside  the  thorax.  Different 
parts  of  the  lung  experience  different  wave  loading  histories  depending  on 
their  geometry,  location  and  neighboring  organs.  At  the  corresponding  location 
where  the  esophageal  pressures  are  taken,  the  I  TP  predictions  have  shown  good 
agreement  with  experimental  measurement.  The  constructed  modal  is  shown  to  be 
an  effective  approach  toward  better  understanding  of  the  blast  biology. 


1.3  qgretAi.  cescupno* 

Section  2  summarizes  the  basics  of  sheep  anatomy.  Section  3  outlines  the 
mathematical  basis  and  the  formulation  of  the  finite  element  method.  Archi¬ 
tecture  of  the  computer  code  used  for  this  program  ia  documented.  Procedures 
for  exchanging  data  between  the  main  finite  eleswnt  analysis  program  and 
graphics  processor  are  also  documented. 

Section  A  describes  a  simplified  three-dimensional  model  defined  by  16 
hexahedral  elements  representing  internal  organs.  The  muscle  and  rib  cage  are 
represented  by  lumped  parameters  which  are  coupled  to  the  finite  element  equa¬ 
tions.  This  model  is  used  to  reproduce  global  ITP  response  with  qualitative 
agreement  under  low  pressure  (<  10  psi)  loading  range.  This  model  is  also  used 
to  study  the  Influence  of  abdominal  air  on  the  ITP  response.  This  application 
confirms  the  previous  experimental  observation  of  Zuckerman  [14]  and  Clemedson 
at  al.  [2,3]  thr-r.  lung  injuries  are  caused  primarily  by  chest  wall  loading. 
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In  Section  5,  a  two-dimensional  model  is  constructed  with  refined  spatial 
and  temporal  resolutions  In  order  to  study  the  wave  propagation  characteris¬ 
tics  Inside  the  thorax.  Material  parametric  studies  are  carried  out  to  Identi¬ 
fy  the  material  parameters  which  are  important  in  the  body's  I  TP  response  to 
external  blast  wave.  The  model  Is  then  used  to  Illustrate  the  spatial  vari¬ 
ation  of  lung  overpressure  histories  when  a  sheep  la  exposed  to  external  blast 
wave  loadings •  It  suggests  thar  tho  interpretation  of  measured  esophageal 
pressure  should  be  carefully  justified  as  progressive  wave  propagation  could 
result  In  different  response  histories  at  different  parts  of  the  lung  and 
thereby  different  damages. 

Following  this  the  model  Is  validated  with  experimental  results  from 
various  loading  cases,  e.g.,  tour  alngle-peak,  three  double-peak,  and  ten 
cases  from  isolmpulse  studies.  Reasonably  good  agreement  is  obtained  among 
comparisons  made  between  the  model  I  TP  prediction  and  the  experimental  esopha¬ 
geal  pressure  measurement. 
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2.  AWAXOMICAL  8TBUCTIHE  OT  SH2EP 


Sheep  ere  used  as  the  prlaary  experimental  subjects  in  the  BOP  program. 
Various  blast  exposure  experiaents  were  planned  in  order  to  understand  the 
bloaechanlcal  responses  and  the  resulting  Injuries.  Concurrently,  a  mathemati¬ 
cal  model  of  sheep  torso  was  constructed  to  simulate  the  transient  interaction 
of  the  anlaal  and  the  blast  wave.  The  aodel  predictions  of  bloaechanlcal  re*' 
sponses  can  be  validated  by  the  experlaental  measurements,  e.g.,  esophageal, 
intrathoracic  pressures,  chest  wall  displacement,  acceleration  histories.  The 
combination  of  anlaal  experiment  and  aodel  simulation  should  identify  the 
injury  criterion  and  understand  the  Injury  aechanlsa.  The  modeling  technique 
can  then  be  used  to  extend  the  results  to  human  and  various  loading  conditions 
for  the  definition  of  a  generalised  safety  guideline  of  blast  exposure. 

To  aodel  the  response  of  a  sheep  torso  under  blast  wave  incidence,  we 
first  need  to  know  its  anatomical  structure.  We  need  to  know  the  shape  and  the 
relative  orientation  of  various  organs  to  define  the  aodel  geometry.  We  also 
need  to  know  their  basic  aechanlcal  properties  to  prescribe  the  structural 
rigidity  and  to  determine  the  resulting  aechanlcal  stress  under  external 
loadings.  Proper  boundary  conditions  should  be  prescribed  to  represent  the 
blast  environment  the  animal  encounters.  With  all  this  information  together 
with  appropriate  engineering  approximations,  a  finite  element  aodel  can  be 
constructed  for  the  animal  torso  to  predict  its  transient  mechanical  inter¬ 
action  with  the  blast  overpressure  wave  loadings. 

2.1  BASIC  SWEEP  AMATOMT 

The  details  of  the  sheep  torso  anatomy  are  summer i red  in  the  following 
illustrations  reproduced  from  References  7  and  9.  Figure  2-1  is  a  side  view  of 
a  sheep  with  superficial  muscle  layers  while  Figure  2-2  is  a  side  view  of  a 
sheep  skeleton.  T^,  L2»  ®tc.  are  labeled  on  the  figures  to  identify  the  loca¬ 
tions  where  a  cross  sectional  view  is  provided.  T^  and  l^,  for  example,  mean 
the  fourth  thoracic  vertebra  and  the  second  lumbar  vertebra.  Cross  sectional 
views  at  T^,  T7,  and  Tg  (Figs.  2-3,  2-4,  and  2-5,  respectively)  show  the  sire 
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15.  ax  oSfiyam  aMOvom  aSdna mo  —  internal  eb- 
Itqua  abdominal  muck 

I T.  at.  ir nolle  rearmA  di  tore  -  thoracic  vontrol 
■emiiH  mutek  * 

IS  ax  aSSyma  name  oSdaaiab  -  eatemal  ob¬ 
lique  abdominal  muck 

IS  qpummat  ad  adPyai  aSdaoioa  fatma  re  tai  ab 
flym  atdaaka  iobnii  —  aponntrotit  af  exter¬ 
nal  and  internal  oblique  abdominal  muckt 

20.  ax  pttaaalk  pi  /nude  fax  pienaoA  atmabmi  - 
deep  pectoral  muck  (aecendinf  pectoral 
mutek) 

21.  ax  Oatair  jhariar  anodmrPa  —  tenter  fatciac  an- 
tebracbli  mutek 


22.  Ityaf  Sayan  ati  fnrlprfa  SmrSit  -  leaf  head  of 
tricepe  bnchii  mutek 

"I.  mpat  bomb  ati  njnjptta  SrarSti  —  tateral  bead 
of  tricepe  braebii  mutek 

24.  ax  SiacSbA  -  brachial  mutclc 

23.  ax  nmat  rorpr  mdaA  -  eatentar  carpi  radl- 
alit  mutek 

26.  ax  r.iOtm  r  mrpt  abara  -  catentor  carpi  ulna- 
rit  muck 

27.  ndrr  (mar  —  tuber  ceaae 

28  ax  tratar  >urmr  Saar  -  tentar  fatoac  latac 
muck 

29.  ax  *Sama  aadlkt  -  middk  (luteal  mutclc 

30'.  ax  pbrndarpi  (30.  ax  (beta  uyr/irta/u,  30'. 
ax  Sirrpi  jhaam>  -  (luteobkept  mutek  (30. 
tupcrficial  (luteal  mutek,  30'.  bkept  fento- 
rit  muck) 

31.  ax  trmtbadmaat  -  temitendinotut  mutclc 

32.  ax  prntarta  Sayo  -  peroneut  lortfut  mutclc 

33.  aSrr  -  udder 


Figure  2-1.  Side  of  superficial  wade  layar  of  a  ahaap 


1.  mmxiUn  —  m»ilN 

2.  miittJikmln  —  mandible 

3.  mtini  —  atlas 

4.  mxk  -  axis 

5.  rertrim  rvnvWu  K  —  fifth  cervical  vertebra 
6-  ter#4m  <k(iciri  /  —  fint  thoracic  vertebra 

7.  ur»6rw  tktvcirn  VII  —  seventh  thoracic  verte¬ 
bra 

3.  wi*h»  /n«4n/a  /  — •  first  lumbar  vertebra 

9.  mtrir*  famhdit  VI  —  tinh  lumbar  vertebra 

10.  mwtrmm  —  uctvm 

1 1 .  ivrirbme  mri^mr  —  coccygeal  vertebrae 


12.  ruM  1 7/  —  seventh  rib 

13.  cmkiXlIt  —  thirteenth  rib 

14.  «rr«f  emtmtit  —  costal  arch 

1 5.  m*/M  tfrmf  —  body  of  sternum 

16.  irqM»  -  scapula 

1 7.  iiMT/w  —  humerus 
li  mhm  —  ulna 

19.  mAm  -  radius 

20.  mm  mrpf  —  carpal  bones 

21.  m  mHmtnrfmlr  III  H  tV  —  third  and  fourth  me¬ 
tacarpal  bone 


22.  m*  dip*****  ***** «  -  phalanges  of  thouc* 
appendage 

23.  m  r*%«*r  —  os  coxae 

24.  m  frmm*  —  femur 

26.  /mtrMn  —  patella 

26.  fifcff  -  tibia 

27.  »a«r  Oai  —  tarsal  hones 

26.  m  mMmton *k  ///r#  /I'  -  thiid  and  fouiOi  n*« 
tarsal  bone 

29.  ttigibMmm  fm*4n  —  phalanges  of  lK  *vkC  * 
pendage 


Figure  2-2 .  Side  view  of  a  aheep  Skeleton 


and  relative  location  of  aajor  organa  in  the  thoracic  cavity*  At  Ty  and  Tg  the 
abdoainal  contents  occupy  aost  of  space.  At  T12  (Fig.  2-6)  we  can  see  a  snail 
portion  of  the  basal  tip  of  the  right  lung.  This  gives  us  an  Idea  of  the  doaed 
shape  of  the  diaphraga  which  contacts  the  lower  boundary  of  the  lungs.  Cross 
sectional  views  at  and  I5  (Figs.  2-7  and  2-8)  show  various  abdoainal  con¬ 
tents  at  different  levels. 

Figure  2-9(a)  shows  the  ventral  view  of  an  isolated  sheep  lung  and  heart. 
Enclosed  in  the  ipsilateral  pleural  sac  each  lung  is  covered  by  the  pulaonary 
pleura  and  is  free  to  aove  in  the  sac.  The  right  lung  [Fig.  2-9(b)]  is  sub¬ 
divided  into  four  lobes  by  interlobar  fissures,  naaely  the  apical t  the  middle, 
the  diaphragmatic,  and  the  accessory.  The  left  lung  [Fig.  2-9(c)l  is  divided 
by  an  Interlobar  fissure  into  two  lobes,  the  apical  and  the  dlaphragaatlc. 
Between  the  lungs  is  the  aediastinua,  soft  tissue  encasing  the  esophagus, 
trachea  and  aorta.  Froa  the  cross  sectional  view  at  and  Ty  we  can  see  the 
esophagus  close  to  the  lung  lobes.  This  is  where  the  experlaental  esophageal 
pressure  is  taken. 

Figure  2-10  shows  the  abdoainal  contents  at  superficial  level  looking 
froa  the  right.  Alaost  three-fourths  of  the  abdoainal  cavity  is  occupied  by 
the  stoaach.  The  stoaach  consists  of  four  coapartaents :  the  ruaen,  reticulum, 
omasum,  and  aboaasa.  The  remainder  of  the  abdominal  cavity  Is  occupied  by  the 
spleen,  liver,  small  intestine  and  large  intestine. 

The  size  of  the  stoaach  and  G1  tracts  is  significant  in  the  sheep  torso* 
It  is  known  to  trap  large  aaounts  of  air  under  normal  circumstances  when 
conpared  with  other  species.  Due  to  differences  in  compressibility  between  air 
and  the  surrounding  tissues,  the  air  has  been  suspected  to  be  directly  related 
to  the  cause  of  GI  tissue  damage  after  animal  blast  exposure.  It  has  also  been 
suspected  of  enhancing  the  possibility  of  lung  injury  due  to  its  presence  in 
the  abdoaen.  A  three-dimensional  aodel  is  constructed  to  study  this  coupling 
effect  in  Section  4. 


* ff*  +  *•  t*  r 
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I.  %  mpmpimb H %.  —ribr  «»  myuyiinl  wd 

MCtwl  UfMVMfMt 

1  ft wot  M—  y»<M  vtHrimr  tkmmrin m  W  —  |pi- 
toui  fNtctu  «f  tixiK  riwfKk  vtntbrt 
1  m  <wf>iwi  (|kM  timmtkm)  -  trtptiiut  immcIc 
•  *h#>r»ck  pan) 

:.J»CW 

1  m.  qpawfc  W  awwppwwife  uti wifii  H  frfrwtl  —  tpi* 
uiU  tt  trmitpiaalU  itoncii  «t  ttfvkii 
mwtcW 

1  m.  ii  yii  II I  m  -  bafiuimw  iktncii 


I 


17.  m  nmmntm  «t  I—  —  iran»«*fw  •!>*»'“ 

fUll  Mtildt 

IB  ml  pri«iill<i  pin  /.  .in  (m  /»f*  uA  Mn4 

inf  protml  punch  pcciuml 

PUMCh) 

It.  m.  puna  m turn.  Mura  UMipfUbina.  r.  fO.pu  w* 
w  -  |Ml  niu.r  Muach.  iifkfuknK  iiml, 
MimyfMurtp 

30l  m,  Ara •  A  nurra  —  urn.  iMiwk  *M 

SI.  V  iiirAnn'imfr  mM  -  cmW  Mimlml 

IS.  piiimi  piitpr  —  Mi  Mu 
11  •  rifbi  Im| 

S4.  MrAMPurairaAb  —  cuubrt  n»rt»wu»M'» 
ti  Mura  rtguib  AaMB  u  fli|ui  —  Bprul  ;ul 
■rub,  nopbpfw 

Ml  M  UhM  |MMu  ffcpA^MUM  —  Ml  —•• 
4UI  cm  *f  Bhpbra(in 

S>.  m.M  i  .Aum  AfMfMI  —  <*mral  •«•• 
<m  *f  ttopbragpi 

SB.  pun  nuMk  AppHi^mulp  -  mill  pun  •*  **' 
‘JB  umu  . . uuurPnA  —  Mnwn  •/  nuuimmi- 


I.  n.  fniuBf*. i  pprAl  fliufifUmlil  -  «pi- 
■ah  «i  «cnhplu»li«  ibnach  «•  rccvici* 

s.  pup  Aram  Mi  ftppr.-N  -  ikmdc  pun  •/  Ira- 
pctiui  Punch 

1  m.  tmftJmm  Burt  -  rtugiiilmm  tbpnwic 
uumii 

4.  Ml  mJyiAu  Miiiuni  —  MuInfLdy!  ibpracit 
until 

V  MM  >1U  Ml  liulin  mm  -  (iftiib  rib,  tr.aier 
punch  W  rib 

B  nulu  Arawu  17//  -  rigbik  ibaruck  v»m  ■ 

bra 

7.  ml  4  n.Hiii  Arab  —  iboruck  Uieratal 
punch 


jy  ^Ml  __  mImmL  g/L  ImMPUMU 

fMM  M,  MM  M7WMM7  •  PMlh  *  M,  NNBTWt* 

t»l  WMielll 

t-  iwNifc^  -  cftftttafe  •(  >c»piil> 

10.  mb  lotoimw  Aw'  —  Intuitu  **r%i  wmirtc 

S3. 

It  ml  tfMfH  miMdlji  Mmiti  —  iktrucic  wuralii 
mraiM  wmc*« 

1 1  Mfe  Xt  —  tWvtwtii  rib 
I  4.  mb  ^Aftm  ftfrnm  mAAmmmmi  •  tMIflHil 
«Miyjt  abOotiwiil  mutcfe 

IS.  MBMB  MfrlWM— K  wi*!*1  *  lift  Xt  — 
ml  mtrtckt.  eleventh  comI  e»rtiUf« 

Id  a.  Hr.  drmmkm  iwfrnMt,  m  mv*m  mSJmmmm  — 
imtrwi  ilwtck  »nm  m0  win,  wraighi  «!»• 
OmuIamI  MMd« 


•It  Bbrf  •  fttitvlif  twlrvi 

32.  «m*mm»  ^brnli  BMMBBrmM  —  rettcvlo -«»•  '*** 
Htcnim 

33.  tMiioPM  MHMLimbMMMMrMMI  —  MMtO.il  •**! 

34.  nwfii»dy»  «  owOmI  n»»  wa 
3S  mSmMMMMMI  ra 

36.  t.  krfmttfm  —  ln^iiic  vet* 

37.  b|Mr  —  liw 

35.  mtmm  n  pnMnr.  Anw  Ar^nnnw  —  h«;unb  •* 
pecirl  vciit,  hepatic  tha 

.1*1.  iiMrM»MM  mtmm  —  Iturr  iMcnOiM 

40.  mimjHh  —  goHliliOOu 

II.  nw»hp  tpbwMbM  —  gtpheid  ctMttttgc 


Vigor*  2-5.  Crooo-ooctlaool  vlw  of  •  ohoop  «t  tho  oighth  thoracic  rortobro. 


II 


1*  * 


Vigor*  2-7.  Cro**-**ctlonal  ala*  of  a  ahaap  at  tha  aacood  lumbar  aartabra 


I.  m.  htftmmm  Mmi  —  iMiflutmu  lamba 


mtm  until 

t  m  nrtybda  t . t  win  —  mnlitfidwr  hiiabn- 


lai-C* 

1  m.  fmmi  muff  -  pilM  najar  r  H»*ck 
l  m.  fmm  km  -  puu  mk».  auucte 
T.  » Mr  Mara tea  -  Imaraal  UUc  anary  and 
aain 

t  a  Mr  tkm mmmt,  »■*' t  -  «mm— I I  Utt 
imr  and  nnln»  ri|hi  irmar 

»  fc.  a— . .  -  mlddla  Utot  lyrapli  aada 

It  m  aMyan  lama  aldiiiaat  -  dwaraal  #4- 

M.  a  kw  jhM  Jam  -  taaaar  faadaa  la«a* 
"Mdi 

It  it  a.  jtin/ifkm id  -  yadrKapr  tia»- 

wamdi 

It  at  aww/ian*  —  wralgti  hmwl  awada 
It  anadwtnmd*  -  laaral  »n«  nmt* 

It.  an  laaka  aBWndha  -  lawnaadlan  «ww 
■nmda 

It  an  laaaa  aM  -  andlal  aawwa  — ada 
•t  a/ mi  ia  a  (taw 


IT.  miQ*a  tana*  aijM  -  laaral  candy* 
af  fmw 

It  trmlftm  aadfcdb  a*  ytaaad  —  kadial  wady* 
af  faamr 

It  an  twyy a^kaanC  —  Uctyi  fawnric  aarit 

Ml  Mria-dbia 

II.  aaaM  —  annilri 

It  an  araa  atdaaM  —  Mralffci  ibdamiail 

M  at  aMfaia  aaiwa  aUani  —  aataraal  afcll- 

H  a  aittifUir  —  addaciar  nyick 
It  ^fw  i  f  J* (  ‘  an  Mia  ayint  —  wyrrfkiil 
laytr  af  |iaur  lamina 


K  (WM ■  pn^mln  aaaad  ai^m*  -  *"»  l»*«  af 
(raaiar  lawwa 
IT.  atM  -  uddar 
It  p/umm  *"  Ja^iaaw 
It  arena  —  cacwca 

30.  mm  tfimSi  ad I  —  tftral  aara  af  calaa 
3|.  artaa  —  nenun 

U  tame  tma  naMaaii  man'  a  —  candadaml 
Wind  rae  af  rvmaa 
U  an  faw*  —  (racilU  armdc 
34.  an  mmtim  -  rarcartui  mutek 
3t  tmfm  a^yaaaa  —  adipaar  bady 
3t  nanruaM  -  kh  wmr 


Vigor*  2-8.  Croat-aact loot  1  vlw  of  a  ihnp  at  tha  fifth  luabar  aartabra 


1.  pau  nuniafai  Mi  ra  liti  -  cranial  part  a f 
cranial  labc 

t  tmmtUit  Mi  Htm'affa  -  caudal  pan  af  cra¬ 
nial  Ufa* 

X  Macmnifan  —  middle  lade 
4.  Mia  rauMa  —  caudal  Mr 
)■  hia erMada  ™  acewierfleie 
&  M|>  dina*l  Mfam3F  -  dactal  margin  Ut- 
Mae)  \ 

7.  may  arena  rra aaffa  -  ventral  mil  margin 
X  /fauna  mwdifaiia  rtaiiiaU  —  cranial  iatetUhar 
fissure 

X/faauafa  awMere  laadMfa  —  caudal  laurlda 
(laaura 

IX  iiiprinr  mmftitm  -  aartic  hripcettion 
II.  hrp*Ji> intdtcrc  -  cardiac  impmtiea 
II  necfaca  —  iracfaca 

IX  a.  pmbmmm/k  drma  —  rlgfct  pulmanacy  artery 

IX  rr.  pul . ala  —  putmcnary  ueiat 

II  Mcrcm  pfaaaar  auCMru m*  -  lawniaa  af  me¬ 
diastinal  pleura 

IX  udhan  carer  waddle  —  ealcul  far  caudal  «ca* 


17.  /faciei  dydayn*a  -  diaphragmatic  cacfacc 
IX  Mntdm  wacdralii  —  tracheal  branches 


(b)  light  long  (fro*  t*f.  2,  Vol.  2, 
fig.  68)a 


I.  pea naaidfa Mi rmuMfa  •.cranial  panaf  CM- 

H|||  ||g| 

X  pan  raadyfa  Mi  rcanaaifa  -  caudal  pact  af  cm- 
nialUfcs 

X  Mm  caadUb —  caudal  labc 
X  may  duet  Miami  -  da  real  margin  Ufa 
mm) 

X  may  arena  rmradfa  —  ventral  acute  marfia 
X  /faaaiu  jam*  tacit  flan  AM  -  imrrtubar  intact 
(caudal) 

7.  may  arena  Main  —  Paul  acute  margin 
.X  faafanecaidtata  -  caediac  natefa 
(  X  Jmrm  dy»ny  mlim  —  diaphragmatic  turf  ace 
rX  IwpcmiwmuCiect  —  aartic  impreuian 
II.  Mipuraim  tupMpa  —  caapfaageal  impeetuna 
IX  iwfih#  wrdkkiiirAi  —  llm  W  m* 
dlexd'ul  pleura 

13.  prmtifmt*  tiukkr  •  left  main  brvn* 

ctma 

14  rr.  pmhmmmSn  -  galiMnary  v*»n» 

11  a  /hJmmmtdrn  tmithm  -  left  pulmonary  rnwy 

II.  frepimei  cwfaw  —  cardiac  impmtion 


(e)  Left  lung  (froa  Kef •  2,  Vol.  2t 
Fig.  67). 


rigure  2-9.  (Cont’d) 
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»  mm  ■actumraim  am  mrmjc  mm  mi  ajcaranal  IgUtH  *»■  .*  utt  an  trwtrnirm  STt  ■  rj,  !"J<fUtXKTi. 


n 


1.  fmki  XU  —  twelfth  Hb 

2.  mm.  mlrnwimin  —  intercostal  muscles 
1  mdhrtrr  «  ri|ht  kidiuy 

4.  pan  tmk  hi  JiejJtrmgmmtk  —  combI  part  of  dia¬ 
phragm 

&  Arpur  —  liver 

5-  tmkm  dnak  —  descending  colon 

7.  pan  mrrmkm  dadw  -  ascending  part  of 
duodenum 


8.  pan  drsrrwdwi  Aidin'  —  descending  part  of 
duodenum 

9.  pan  rmaMij  dasdmt  mira  JHIrm  —  cranial  part 
of  duodenum,  gallbladder 

10.  rwrtjlmgm  tmmki  X  -  tenth  costal  cartilage 

11.  arrtrr  dkrftrr  —  right  ureter 

12.  mmmipimtn  rW»  —  spiral  ansa  of  colon 
El  nr  aw  —  cecum 

1 4.  amapntrMudti  raft  -  proaimal  ansa  of  colon 


1 5.  p/nh«m  —  jejunum 

*6.  parwv  mfirrfirinlu  tmrmti  mm/mi  —  superficial 
layer  of  greater  omentum 

1 7.  para*  pntfttmlm  mrnmtii  tmjnrit  —  deep  la>’*r  ° 
greater  omentum 

1 8.  MMTN/MN  MMM  —  IctSCr  omCIttUIT 

19.  bImmimw  —  abomasum 

20.  pan rim  —  pancreas 


figure  2-10.  The  abdominal  content*  at  superficial  level 
looking  froa  the  right. 


| 

I 
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i  I'ULrrMHtJUtA  ruinc  ’  ur  1 


L/w^AftAftwu»  r*,'wusA',.\n/i^JWW\i  i*. 


3.1  HAnmmcAL  basis  abs  the  fomulatmm  or  finite  element  method 


3.1.1  Introdoctlon 

Finite  element  method  (FEN)  Is  used  to  construct  a  mathematical  model  ot 
the  animal  body  and  to  study  its  structural  response  under  external  blast 
overpressure.  The  idea  of  FEM  is  to  divide  a  continuous  body  into  a  number  of 
elements  of  finite  size.  The  original  continuous  body  is  replaced  by  t  k-  | 

assembly  of  these  elements.  The  continuous  physical  quantities  are  then  E 

approximated  by  functions  which  are  smooth  in  each  element,  but  are  continuous 
and  only  piecewise  smooth  in  the  whole  body  assembly.  The  governing  differen¬ 
tial  equation  is  solved  in  a  "weak  form"  in  that  it  is  satisfied  in  a  weighted 
average  sense.  The  choice  of  approximations  is,  therefore,  a  key  to  the  FEM. 

Basically,  we  approximate  the  physical  quantities  in  terms  of  nodal  unknown 
parameters  through  spatial  Interpolation  function.  Instead  of  the  original 
differential  equations  we  then  have  a  set  of  algebraic  equations  in  general¬ 
ized  coordinates  defined  at  nodal  points.  The  distribution  in  an  element  is 
obtained  through  interpolation. 

There  is  a  wide  variety  of  approaches  to  FEM  discretization  depending  on 
the  variational  statement  or  weight  residual  method  [13].  Here  a  compatible 
displacement  model  based  on  the  principle  of  virtual  work  is  used.  The  princi¬ 
ple  of  virtual  work  states  that  for  a  body  of  volume  V  bounded  by  surface  A 
the  virtual  work  done  by  the  external  body  force  F^,  inertial  force  and 

surface  traction  force  can  be  written  at  any  instant  of  time  as: 

/  o^d*^  dV  -  /  (Ft  -  d V  +  /  6.^  dA  ,  (1) 

V  V  A^ 

where  and  denote  the  Cauchy  stress  and  small  strain  tensors  and  A^ 
denotes  the  portion  of  the  boundary  surface  where  traction  is  prescribed. 

The  immediate  step  is  to  find  displacement  function  u^  such  that 
on  Au  (displacement  or  Dirichlet  type  boundary  conditions)  and  Eq.  (1)  is  sat-  j 

lsfled  for  all  admissible  virtual  displacement  6uj.  If  the  choice  of  u^  is 
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sufficiently  smooth  the  solution  to  Eq.  (1)  could  well  be  exactly  that  to  its 
original  differential  forn  of  equation  of  motion.  One  important  merit  with  FEM 
lies  in  that  it  requires  the  choice  of  only  to  make  Eq.  (1)  integrable.  In 
other  words,  requirements  on  the  choice  of  admissible  function  is  not  so 
rigid  as  to  satisfy  the  original  differential  equation  exactly  but  only  need 
to  satisfy  it  in  a  sense  of  weighted  average.  This  way  we  gain  flexibility  in 
coping  with  irregular  geometry  and  boundary  conditions.  A  sparse  and  symmetri¬ 
cally  banded  matrix  equation  is  formed  to  be  solved  for  nodal  displacement  at 
each  time  step.  The  body  configuration  is  updated  for  the  Incremental  solution 
at  each  time  increment  before  marching  to  the  next  time  step. 

3.1*2  Spatial  Discretisation 

The  finite  element  discretization  is  done  by  first  dividing  the  whole 
body  into  a  finite  number  of  sub-domains,  each  called  an  element.  For  a  gen¬ 
eral  continuum  model  .todal  displacements  are  chosen  to  be  the  primary  unknown 
to  solve  for.  For  a  3D  model  each  node  is  assigned  with  three  degrees  of 
freedom.  By  looking  at  the  variational  statement,  Eq.  (1),  a  Cq  continuous 
type  interpolation  function  will  be  sufficient  to  represent  the  displacement 
field  in  each  element.  We,  therefore,  write 

ND 

Uj^.y**)  ■  l  Na(x,y,z)  U®  (2) 
a-1 

for  each  element  where  N*(x,y,z)  is  the  element  interpolation  function  associ¬ 
ated  with  node  a  and  is  the  lt^1  component  of  nodal  displacement  at  node  a. 
Note  the  summation  in  Eq.  (2)  is  to  be  carried  up  to  the  total  nodal  number  ND 
associated  with  an  element. 

In  our  application  we  shall  construct  both  3D  and  2D  models  for  different 
purposes.  In  the  3D  model,  eight  node  brick  elements  are  used  with  trillnear 
interpolation  functions.  In  the  2D  model,  4-node  quadrilateral  elements  are 
used  with  bilinear  interpolation  functions.  Isoparametric  formulation  is  used 
in  the  formulations: 


With  substitution  of  Eqs.  (2)  and  (3)  into  Eq.  (1),  we  yield: 


ab 


*  ("tJ 

e-1 


Ui  + 


*5  -  $  - 


(4) 


where 


-  /  P  N*Nb  dV  6.. 

1J  v*  ° 

is  the  consistent  mass  matrix; 


/. 

V* 


1J 


is  the  Internal  force  vector;  and 


(5) 


(6) 


R t  ¥  Nb  d*  +  /  p  F,  Nb  dV  (7) 

J  J  ve  °  3 

•• 

is  the  generalised  force  vector.  The  first  term  in  Eq.  (4),  with  denoting 
nodal  acceleration,  is  the  inertial  force.  The  second  term,  the  internal  force 
vector,  accounting  for  the  force  due  to  elasticity,  can  further  be  written  as 
the  sum  of  contributions  from  Incremental  stress  and  initial  stress.  The  last 
term,  or  generalised  force  vector,  specifies  the  generalised  nodal  force  due 
to  prescribed  body  force  and  surface  tractions.  Note  that  Kb,  and  Rb  are 

the  corresponding  matrices  or  vector  in  element  level.  The  integration  is  car¬ 
ried  at  element  level  with  V®  and  Ae  denoting  element  volume  and  surface  area, 
respectively.  The  global  form  of  dlscretlsed  equation  of  motion  is  obtained 
through  assembly  of  contributions  from  all  elements.  The  assembly  procedure  is 
symbolically  written  in  terms  of  summation  through  total  number  of  element  E. 

At  a  given  time,  once  the  global  equation  of  motion  [Eq.  (4)]  is  assem¬ 
bled,  it  is  then  modified  for  non-sero  Dirichlet  boundary  conditions  and/or 
general  mixed  type  boundary  conditions.  Note  that  Neumann  type  boundary  condi¬ 
tion  enters  the  problem  by  dictating  the  variational  statement,  Eq.  (1). 


3.1.3  taporri  Mgcretl««tloa 


We  next  have  to  do  temporal  discretization  of  Eq.  (4)  since  It  Involves 
time  differentiation.  A  direct  integration  scheme  based  on  Newmark  approxima¬ 
tion  Is  used: 

Vl  ■  °n  +  “  ”n  +  At'  [(?  -  B)°„  +  8  Vl]  «> 

v.  -  K  *  “  [(1  -  *>  K  * '  Vi]  <»> 

In  which  the  subscripts  denote  variables  at  cn+l  and  tn>  respectively,  and  At 

“  cn+l  “  tn*  ^  and  U  are  nodal  velocity  and  acceleration,  respectively.  0  and 

y  are  parameters  that  control  numerical  damping  and  stability.  For  linear 
problem  y  ■  1/2  produces  no  numerical  damping.  By  choosing  20  >  y  >  1/2  uncon¬ 
ditional  stability  Is  achieved  regardless  of  the  size  of  time  step  At. 

At  time  tn+1  we  can  rewrite  the  glogal  form  of  generalized  equation  of 
motion  In  matrix  notation  as: 


»  Vi +  W  ■  Vi 


(10) 


Substi  ion  of  Eqs.  (8)  and  (9)  into  (10)  yields 


0At 


2  Vi ♦  <  Vi)  *  *. 


L  +  A 
n+1  n 


-  1  . .  + - y  0  +  At  D 

n+1  0At2  I  n 


2  ”  / 
01  At  Unj 


.  (11) 


To  solve  Eq.  (11)  for  Un+j  we  can  first  rewrite  It  In  a  different  form, 


*Vi)  ■  0  • 


(12) 


The  question  Is  then  to  find  solution  such  that  Eq.  (12)  is  satisfied 

within  specified  error  tolerance.  We  may  write  Taylor  expansion  of  F  about 
«S+1  where  U^+j  denotes  an  Initial  test  point: 
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♦  40n+l)  ■  *0  +  "«».)  •  4D-*.  +  <hl«h  “J«  “"*•>  • 


with  D  denoting  Prechet  derivative  operator*  Iteration  formula  can  be  obtained 
with  linearization  of  Eq.  (13): 


“0>L)  •  iDi*i  -  - 


•it!  ■  °i+i  +  4Di+1  •  <>5> 

For  a  general  nonlinear  form  of  Bq.  (12),  Iteration  la  repeated  until  AU^,,  is 

n+l 

smaller  than  certain  specified  error  tolerance* 

Following  these  steps  Eq.  (11)  Is  to  be  solved  from 


+r-?i  •  i0i+i  -  t*n«  -  -  “111 


•it!  ■  °i+i +  40i+i  <17> 

Equation  (16)  can  be  rewritten  as 

•  iDl+l  '  *L  <lf,) 

A  A 

with  K  and  *.  denoting  effective  stiffness  matrix  and  effective  load  vector, 
respectively* 


3*1*4  Material  Description 

In  our  application  the  materials  of  various  organs  are  assumed  to  be  iso¬ 
tropic,  linear  viscoelastic.  The  constitutive  relationship  is  written  as: 

olj(t)  -  J*  6  [«t  -  t)  -  |  G(t  -  t)]  de^  +  2  /  C(t  -  t)  de^  (19) 
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where  Ojj  and  denote  Cauchy  etrees  end  small  strain  tensors ,  1,  G  are  bulk 
and  shear  relaxation  moduli,  respectively,  and  6^^  Is  the  Kronecker  delta 
function.  If  we  further  assume  that  the  body  la  completely  undisturbed  prior 
to  time  0  we  can  rewrite  Eq.  (19)  as 

°lj(t)  *  +  /\j  [K(t  -  O  -  f  G(t  -  T)j  dEkk  +  2  J* C(t  -  t)  dctj  (20) 

where 

°lj(t)  ’  6lj  [K(t>  "  I  ekk(o)  +  2C(t)  elj(o)  * 


For  practical  application  the  relaxation  moduli  are  approximated  by  the 
following  exponential  series  (Ref.  11): 


m  -t/x 
K(t)  -  Ko  +  l  exp 


(21) 


and 


K  -t/B 

G(t)  -  G  ♦  l  G.  exp  1  . 


L 

J-l 3 

Substitution  of  Eqs.  (21)  and  (22)  Into  (20)  gives 


-t/X  t  t/X 

d  r  n  , 

*  '  -y  A‘"  *  I  6 


o1:J(t) » oj4  +  «44  i  *. 


(22) 


kk 


1C  -t/B_  t  t/B  /  2  \ 

+  J^n  e  l  e  (M«1J  -T6lj  dekk)  ’  <23> 

If  we  assume  the  strain  varies  linearly  within  each  time  increment  we  can 
write 


dt 
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and  eha  resulting  Integrals  can  be  evaluated  to  give  a  recursion  formula  for 
each  term  In  the  sum.  Two  forma  are  used: 


-At/X  X  -At/X 

*VCJ  ♦'.ife-*  *) 


(24*) 


and 


Vv«> 


-At/e. 


P_  — At/B  r  o  l 

+  Cp  At  ( 1  "  *  )  [2Aclj  “  T  6ijAekk]  (24b) 


where 


it  »  t  . .  -  t 
n+1  n 


and 


Aekk  "  ekl^tn+l^  “  ekk^n^  * 

In  terms  of  and  Jp  we  can  rewrite  Eq.  (23)  as: 

°lj^tn+!^  "  6lj  (*o  “  1  Co)Ekk^tn+l^  +  2Coelj^n+l^ 


+  6.,  i  I(txl)  +  V  J  ft  A 
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The  tangent  modulus  at  tn4-1  is  written  as: 


17  ,  \  *  X  -At/X  -  K  e  "At/Bl 

“IJM  ■  *ij*u  [K  -  T  eo)  +  !?('-•  )  •  f  i  GP  if 
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This  Is  the  expression  of  DK  used  in  Eq.  (16). 


3.2  OOKFUTKl  COM  AKCH1TVCTTJKE 

Based  on  the  above  framework  the  FEM  is  used  as  a  tool  to  study  the 
transient  response  of  a  sheep  torso  model  under  external  blast  overpressure. 
Consideration  of  sheep  anatomy  is  described  in  Section  2.  Materials  of  various 
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organa  art  assumed  to  ba  doacrlbad  by  tha  laotroplc  linatr  viaeoalaatlc  law  aa 
out 1 load  In  Sactlon  3.1.4. 

Tha  computer  coda  FEAPWS  la  uaad  to  parfora  thla  analyala.  It  la  a  var- 
alon  modified  froa  FKAP74  developed  by  1.  L.  Taylor  of  UC  Barkalay  [11).  It 
waa  adaptad  to  run  on  tha  JAYCOB  aalnfraaa  Uni vac  1100/60.  Tha  atoraga 
arrangaaant  haa  baan  raatructurad  for  clarity  and  uaar  control.  Tha  prograa 
haa  baan  aodlflad  for  dlffarant  apadflc  naada  and  taaka. 

Tha  computer  coda  FBAFWS  conalata  of  two  functional  units:  Data  Input 
nodula  and  aquation  noising  and  output  aodula.  In  data  Input  nodula  tha  uaar 
daflnaa  hla  apadflc  boundary-inl tial-valua  problem  by  supplying  tha  following 
Information  (Fig.  3-la): 

1.  Control  parameters  -  Including  Input  of 

Dimensions  -  slsa  of  tha  problem 

Humber  of  degrees  of  freedom  par  node 

Total  number  of  elements,  nodal  points  and  materials 

Dimensioning  slsa  of  storage  arrays 

Total  number  of  time  staps,  time  Increment  of  each  time  step 
Numerical  method,  stability  and  damping  control  parameters 
Output  data  managing  parameters 
etc. 

2.  Material  data  -  input  of 

Mechanical  properties 

Physical  properties 

Kind  of  structural  element  to  use 

etc. 

3.  Nodal  data  -  for  each  node.  Input  of 

Nodal  number 
Coordinates 

Boundary  constraint  code 
etc. 

4.  Element  data  -  for  each  element.  Input  of 

Element  number 
Element  connection  array 
Material  type 
etc* 

3.  Boundary  conditions  -  specification  of  magnitude  and  direction  of 
displacement  boundary  conditions  In  terms  of  generalised  nodal  dis¬ 
placement  . 

6.  Force  loadings  -  specifications  of  magnitude  and  direction  of  point 
nodal  force,  surface  traction  and  body  force.  These  loading  are  all 
converted  Into  equivalent  generalised  nodal  force. 
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Figure  3-1 (a) .  Flow  chart  showing  data  Input  nodule 


Djmaalc  irinilant  analyata 
(taaUctt  lata* ration) 


showing  solution 
module . 


o~ 


t  •  t  +  At 


wr  •  nr  +  i 


Figure  3“l(b).  (Cont'd) 
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4cij  -{(‘"l.J  +4“j.t) 

°ij  “  f(Eij’  '  (Eq*  25) 

(B)  Output  strain  and  stress  components,  e^j(t)  and  o^j(t). 

(C)  Compute  the  effective  load  vector  Rj(t  +  At)  (R.R.S.  of  Eq.  18). 

(D)  If  NT  *  1  (not  the  first  time  step)  go  to  (E). 

(a)  Form  diagonal  lumped  mass  matrix 


(b)  Fora  effective  element  stiffness  matrix 

I*  -  M6  +  EK  . 

SAtZ 

(c)  Apply  displacement  boundary  condition* 


(L.H.S.  of  3q.  18) 


le 

(d)  Store  entries  of  a  in  global  effective  stiffness  matrix  K 

(Fortran  name  A(I,J)) 

Ae  A 

(e)  Store  entries  of  R  in  global  effective  load  vector  R. 

(Fortran  name  DF) 

*  m 

(f)  Decompose  K  »  LDL 

(g)  Go  to  4 

(E)  (a)  Apply  displacement  boundary  conditions,  if  needed* 

A  g  A 

(b)  Store  entries  of  R  in  global  effective  load  vector  R* 

(Fortran  name  DF) 


(c)  Retrieval  of  L  and  D. 

A  A 

4.  Solve  Au  from  K(t  +  At)  •  Au(t  +  At)  -  R(t  +  At) 

5.  Update 

u^(t  +  At)  -  u^(t)  +  Au^(t  +  At) 


(Eq.  18) 


(Eq.  17) 


u4(t  +  At)  -  (l  -  j)  u±(t)  +|l  -  -|^At  Ui(t)  +  (gk)Au1(t  +  At) 

(Eq*  8) 

Ui(t  +  At>  "  -fest)  “i(t)  +(J  "  h)  Ui<0  +  (“l)*1 i(t  + 


(Eq.  9) 


6.  t  -  t  +  At 
NT  -  NT  +  1 

Go  to  1  for  next  time  step. 


For  brevity  these  steps  are  summarized  in  the  flowchart  on  Figure  3-1 (b). 
Note  that  this  algorithm  is  for  linear  problems  only.  For  nonlinear  problems 
iteration  should  be  called  for  to  obtain  Au  under  prescribed  error  tolerance 
from  Eq.  (18)  in  step  4. 
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3.3  DATA  EXCHANGE  BETWEEN  FINITE  ELEMENT  ANALYSIS  PROGRAM  AND  GRAPHICS 
PROCESSOR 

MOVIE  is  an  interactive  graphics  program  installed  in  the  standalone 
SUPERSET  graphics  processor  (SST)  while  FEAPWS  is  the  finite  element  analysis 
program  on  the  Uni vac  1100/60.  The  linkage  between  these  two  codes  will  enable 
us 


•  To  define  and  modify  the  geometry  of  finite  element  models  with  MOVIE 
and  at  the  same  time  display  it  on  a  graphics  terminal  until  satis¬ 
fied, 

•  To  visualize  the  body  deformation  and  display  various  function  dis¬ 
tributions  of  the  model  under  analysis. 


To  achieve  this  the  following  steps  must  be  followed  (see  Figure  3-2) : 

1.  Define  the  model  geometry  with  MOVIE  in  SST. 

2.  Rewrite  the  geometry  data  in  the  "Univac  way”  since  the  ASCII  code  is 
written  in  different  ways  on  different  machines.  This  is  done  with 
UTL. RE AD80/ ASCII. 

3.  Once  the  data  ia  transmitted  to  the  Univac,  rewrite  the  geometry  file 
in  FEAP  format.  Information  on  material  properties,  loadings,  bound¬ 
ary  conditions,  initial  conditions  and  the  FEM  model  basic  parameters 
as  well  as  the  numerical  scheme  parameters  are  added,  then  an  input 
file  for  FEM  analysis  is  prepared.  This  subtask  is  performed  by 
LINK.MOVFEA. 

4.  Upon  completion  of  FEM  transient  analysis  with  FEAPWS,  the  computa¬ 
tional  results  (displacement,  stress  components,  etc.)  are  written  to 
an  output  file  in  MOVIE  format.  This  is  done  by  LINK.FEAMOV. 

5.  Rewrite  the  resulting  file  in  "SST  way"  with  UTL.WRITE80/ASCII  so  it 
is  recognizable  to  SST. 

6.  Once  the  data  is  transmitted  to  the  SST,  then  the  graphics  display 
functions  of  MOVIE  can  be  used  to  visualize  the  computational  re¬ 
sults.  Even  then,  the  procedure  to  display  any  result  is  quite  tedi¬ 
ous.  The  results  are  written  to  separate  data  files,  e.g.,  geometry, 
displacement  (nodal  values),  and  scalar  functions  (element  values), 
for  each  time  step.  A  long  sequence  of  Interactive  commands  have  to 
be  supplied  by  the  user  to  manipulate  the  data  into  a  form  ready  for 
graphic  display.  To  simplify  this  procedure,  file  manipulation  com¬ 
mands  and  interactive  display  commands  have  been  packed  and  written 
as  two  programs,  FILEPREP  and  LOOK  which  are  executed  by  the  user  at 
the  SST  terminal. 
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Figure  3-2.  Block  diagram  summarising  atop*  In  linkage  automation 
between  MOVIE  and  FEAFV8. 
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7.  If  further  corrections  or  modif ications  (spatial  element  refinement, 
varying  geometric  size,  shape,  material  properties,  different  load¬ 
ings,  boundary  conditions,  etc*)  are  necessary,  steps  1  through  6  are 
repeated* 

Programs  written  for  each  of  the  specific  steps  in  this  section  are  re¬ 
produced  in  Appendix  A  for  future  reference. 
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4.  A  SIMPLIFIED  THKEE-DIMEMSIOHAL  MODEL 


4.1  BBSCBIPnOE  OV  THE  MODEL 

4.1.1  trMtgwU  of  lap<d  Material  Frepertif  on  the  FPC  Modtl 

Based  on  the  sheep  torso  anatomy,  a  three-dimensional  FEM  model  of  16 
eight-node  brick  elements  and  64  nodal  points  is  defined  to  capture  the  major 
organ  responses  important  in  blast  overpressure  exposure  of  a  sheep.  Figure 
4-l(&)  is  an  exploded  view  shoving  representations  of  various  major  internal 
organs  in  the  thorax  and  abdomen.  Figure  4-1 (b)  indicates  the  dimensions  of 
different  parts  in  the  model. 

To  describe  all  the  other  structural  aspects  of  the  anatomy  (muscle 
layer,  rib  cage,  spinal  column,  etc.)  and  not  Increase  the  number  of  solid 
elements  used  in  the  analysis,  it  is  necessary  to  implement  a  model  to  combine 
lumped  material  properties  (suggested  by  Professor  Y.  C.  Fung).  The  idea  is  to 
lump  the  inertial  and  elastic  properties  of  the  rib  cage  and  muscle  layers  to 
the  edges  and  nodal  points  of  those  elements  describing  the  Internal  organs. 
This  treatment  brings  the  contributions  of  rib  cage  and  muscle  layer  into 
consideration.  It,  however,  greatly  reduces  the  size  of  matrix  equation  to  be 
solved  at  each  time  step  and  makes  the  computational  cost  and  time  feasible. 

To  explain  this  in  detail  we  shall  write  the  global  linearized  form  of 
equation  of  motion  for  internal  organs  as  [Eq.  (16)  of  Sec.  3.1.3] 


where  M  represents  the  contribution  due  to  mass  inertia  of  various  internal 
organs.  The  additions  of  muscle  layer  and  rib  cage  are  justified  as  follows: 

Muscle  Layer  -  The  muscle  layer  is  considered  to  be  "bordering"  the 
external  surface  of  the  internal  organ  medium.  To  lump  its  mass  to  the  FEM 
representation  of  the  internal  organs  we  first  have  to  determine  the  volume  of 
each  of  the  muscle  masses.  This  is  estimated  by  multiplying  the  surface  area, 
obtained  from  the  coordinates  of  the  4  nodal  points,  by  an  average  thickness 


(a)  Exploded  view  showing  various  organs  represented 


Figure  4-1.  Sheep  torso  model  (head  up  posture) 
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of  3.5  cm.  The  density  Is  assumed  to  be  1  gm/ca.  The  total  mass  of  a  muscle 
layer  patch  Is  to  be  shared  by  the  four  nodal  points  defining  the  bordering 
surface.  Therefore,  one-quarter  of  the  mass  of  this  muscle  patch  Is  added  to 
the  corresponding  diagonal  term  of  the  mass  matrix  N  In  Eq.  (16).  In  the 
computer  program,  this  Is  done  at  the  element  level  before  the  assembly  of  the 
global  equation. 

lib  Cage  -  The  estimation  of  lumped  mass  due  to  rib  cage  Is  done  by  first 
considering  that  only  a  portion  of  the  external  surface  of  lung  elements 
borders  the  rib  cage.  This  effective  bordering  surface  Is  hypothetically 
assumed  to  be  one-third  of  the  surface.  In  other  words,  for  a  lung  external 
surface,  one-third  of  the  area  Is  under  the  bony  rib  which  has  a  thickness  of 
1  cm  and  mass  density  of  2.8  gm/cm  •  Its  total  mass  Is  then  added  to  the  four 
nodal  points  defining  the  external  bordering  surface.  Computationally,  one- 
quarter  of  this  mass  Is  added  to  the  corresponding  diagonal  term  of  the  mass 
matrix  N  In  Eq.  (16). 

When  the  body  Is  under  external  surface  loading,  the  stiffness  of  the 
whole  rib  cage  arises  from  the  sMffness  of  the  connecting  ligaments  and  the 
bending  rigidity  of  individual  bony  elements.  To  study  the  global  character¬ 
istic  responses  of  the  lung  under  blast  wave  loadings,  a  detailed  analysis 
accounting  for  the  bending,  compression,  and  shear  deformations  of  the  rib 
cage  does  not  seem  to  be  beneficial  at  this  stage. 

Similar  to  the  treatments  of  lumped  mass  for  the  muscle  layer  and  rib 
cage,  the  structural  rigidity  of  the  rib  cage  Is  lumped  to  the  finite  element 
stiffness  representation  of  the  surface  Internal  organs.  The  procedure  Is  to 
lump  a  structural  frame,  which  consists  of  bars  of  high  stiffness  In  its 
direction  of  axial  compression  and  angular  rotation,  such  that  edges  defining 
the  lung  elements  In  the  sheep  model  move  in  conjunction  with  the  "rib  cage” 
under  any  external  force  loading  (see  Fig.  4-2).  The  stiffness  ot  bars 
constituting  the  rib  cage  Is  chosen  to  be  1.7  x  10  dyne/cm  .  The 
contribution  of  stiffness  due  to  thl?  lumped  rib  cage  Is  added  to  the 
corresponding  force-deflection  constants  and  DK^j  of  Eq.  (16). 
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(a)  Structural  frane  repraaentlng  rib  cat**  Thla  fr«u 
con* lata  of  hitfi  axial  atlffnaaa  bara  and  will  cap 
the  torso . 


(b)  Slow  of  terao  nodal  with  elrelta  denoting  luapod  aaaa 


Figure  4-2.  Lumped  rib  cage  wodel 


Table  4-1*  Material  Propertlaa  of  Varioua  Organa  Uaad  In  tha 
Simplified  3D  Finite  Element  Model 


Organ 

Lung 


Rumen 


Omasum 


Small 

intestine 


Large 

intestine 


Property 


*1 
» 1 
P 


Value 

Source 

1.4  *  10®  dyne/c*2 

Ref.  6 

1.0  *  JO®  dyne/cm2 

computed  from 
t3R(l-2v)]/[2(l+v)] 

0.30  .  , 

Ref.  6 

0.7  x  10®  dyne/cm2 

assumed 

0.5  x  10®  dyne/cm2 

assumed 

1  msec 

assumed 

1  msec. 

assumed 

0.37  g/cm3  (small  animal) 

Ref.  5 

0.19  g/cm 3  (large  animal) 

Ref.  5 

1.0  x  A.O10  dyne/cm2 

1.0  x  10®  dyne/cm2 

assumed 

Ref.  12 

0.50  o  * 

assumed 

5  x  10®  dyne/cm2 

assumed 

5  x  103  dyne/cm2 

assumed 

1  msec 

assumed 

1  msec 

assumed 

1.0  g/cm3 

assumed 

1.0  x  10}°  dyne/cm2 

1.0  x  10®  dyne/cm2 

assumed 

Ref.  12 

0.50 

aaeumed 

5  x  10®  dyne /cm2 

assumed 

5  x  103  dyne/cm2 

assumed 

1  msec 

assumed 

1  msec 

assumed 

1.0  g/cm 3 

assumed 

1.0  x  10*°  dyne/cm2 

2.6  x  105  dyne/cm2 

0.50  Q 

5  x  10®  dyne/cm2 

1.3  x  103  dyne/cm2 

assumed 

Ref.  12 
assumed 
assumed 
assumed 

1  msec 

assumed 

1  msec 

assumed 

■**  t* 0  g/cmJ 

assumed 

1.0  x  ioJ°  dyne/cm2 

7.5  x  10®  dyne/cm2 

assumed 

Ref.  12 

0.50 

assumed 

5  x  10®  dyne/cm2 

assumed 

3.8  x  10®  dyne /cm2 

assumed 

0.5  msec 

assumed 

0.5  msec 

assumed 

1.0  g/cm 3 

assumed 
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4.2*1  81mgls  I»k  loowrt 


Pour  cases  of  alngla  peak  loading  are  demonstrated  in  Figure  4-3.  The 
four  experiments  were  performed  on  different  animals  under  different  amounts 
of  TNT  frea-fleld  explosion  at  various  distances  fro*  ths  explosive  center. 

Bach  of  the  experimental  esophageal  pressure  recordings  shows  the  local 
lung  tissue  Is  under  compression  phase  during  the  first  10  ms  with  peak  pres¬ 
sure  occurring  at  4-7  ms  from  the  onset  of  the  blast  arrival  at  the  model  body 
surface.  Approximately  10  as  after  the  onset  the  local  lung  tissue  starts 
experiencing  tensile  phase  with  slower  response  and  reduced  amplitude.  Reduc¬ 
tion  In  amplitude  and  changes  In  frequency  may  be  due  to  the  damping  occurring 
Inside  the  body. 

The  16-element  body  representation  la  certainly  a  simplified  model.  The 
purpose  Is  to  capture  the  main  body  reaponae  characteristics.  For  that  con¬ 
sideration  ITP  predictions  of  different  cases  (Fig.  4-3)  demonstrate  reason¬ 
ably  good  agreement  with  experimental  measured  esophageal  pressures  when  the 
blast  Is  less  than  10  psl. 

4.2.2  Pomble  Beak  Exposure 

Blast  pressure  recordings  taken  from  the  double  peak  study  of  Lovelace 
ITRI  during  October  1981  and  May  1982  are  also  applied  on  the  3-D  model  to 
study  Its  Intrathoraclc  pressure  responses.  Figure  4-4  shows  comparisons  of 
ITP  prediction  and  experimental  measured  esophageal  pressure  for  three  cases 
with  different  time  duration  between  two  blast  peaks.  These  comparisons  shows 
that: 

1.  For  At  -  3.6  and  7.6  as  the  two  humps  of  the  calculated  ITP  response 
to  double  peak  blast  are  not  distinct  as  are  the  data. 

2.  For  At  ■  9.7  as  the  timing  and  low  frequency  response  of  model  pre¬ 
diction  are  comparable  to  those  of  experimental  measurements.  But  as 
In  the  other  two  cases,  the  wave  fora  of  the  experimental  measured 
esophageal  pressure  has  a  faster  rise  and  larger  peak  pressure. 

For  these  cases  the  peak  blast  rises  to  40  psl  In  0.5  ms;  It  also  falls 
off  to  0  within  2  to  3  ms.  It  appears  that  a  more  refined  model  is  needed  to 
provide  both  spatial  and  temporal  resolution  under  these  ranges  of  loading. 
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IMTRA  THORACIC  PRESSURE  (PSII  IMTRfl  THORACIC  PRESSURE  IPSII 

B  -5.0  0.S  5.0  M.B  15.0  30.0  35.0  -OJ-5J  OlO  5.0  *■  H.0  3010 


SHEEP  MODEL 


SHEEP  MODEL 
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TIME  CHS) 

(a)  8  lb  TNT  (free  field), 

range  21.3  ft,  HOB  3.2  ft 


.0  5.0  10.0  IS.O  20.0  25.0  30.0 

TIME  ins) 

(b)  16  lb  TNT  (free  field), 

range  27.0  ft,  HOB  6.0  ft 


SHEEP  MODEL 


SHEEP  MODEL 
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61  lb  TNT 
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Figure  4-3.  Comparisons  of  FEM  predictions  and  experimental  ITP  under 
different  single-peak  blast  loadings. 
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INTRO  THORACIC  PRESSURE  IPS1) 


<•)  At  -  3.6  m 


(b)  At  -  7.6  m 


(e)  At  -  9.7  m 

Figure  4-4.  Comparisons  of  FEM  predl"':ions  and  experimental  ITP 
under  different  double-peak  blast  loadings.  At  is 
the  time  duration  between  two  blast  peaks. 
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4.3  3WSHTVITT  STUDIES 


4.3.1  Influence  of  Abdominal  Air  on  ITP  Response 

Sheep  are  used  as  experimental  subjects  at  WKAIR  and  Lovelace  ITRI  In 
order  to  correlate  the  external  blast  parameter  with  Injury  map.  Sheep, 
hovrever,  have  large  amounts  of  air  In  their  61  tracts.  The  question  of  how 
much  Influence  this  trapped  air  has  on  ITP  response,  therefore,  needs 
clarification. 

As  an  attempt  to  answer  this  question,  the  simplified  16-element  sheep 
model  with  different  degrees  of  abdominal  air-fllllng.  was  exposed  to  blast 
loading.  Figure  4-5  shows  the  results  of  ITP  response  for  various  cases.  In 
case  (a)  the  model  does  not  have  any  air  in  its  abdomen.  In  cases  (b)-(d)  the 
models  are  given  Increasing  amounts  of  air  In  the  rumen.  Case  (d)  has  half  of 
the  rumen/ small  Intestine  region  filled  with  air*  Compared  with  case  (a),  the 
peak  ITP  Increase  Is  less  than  10Z.  It  indicates  that  the  ITP  response  Is 
relatively  insensitive  to  the  abdominal  air.  The  ITP  response  Is  primarily 
caused  by  loadng  on  the  chest  wall. 

This  supports  the  experimental  work  of  Zuckerman  [14]  and  Clemedson  and 
J*dns8on  [2,3]  in  which  they  demonstrate  with  rabbits,  dogs  and  other  species 
that  lung  injuries  are  caused  primarily  by  blast  loading  on  the  chest  wall.  It 
is  noted  that  this  does  not  exclude  the  possibility  of  increased  61  injury 
occurrence  and  severity  on  sheep  due  to  abdominal  wall  blast  loading  when 
there  is  an  excess  amount  of  air  trapped  in  the  61  tract. 

4.3.2  Iaolspulme  Studies 

To  further  examine  the  behavior  of  the  current  version  of  the  model,  iso- 
impulse  studies  have  been  done  in  the  following  two  ways  (rabbit  size  model): 

1.  Isolmpulse  blast  of  10  psi-ms  with  different  dP/dt  is  applied  on  the 
surface  of  the  model.  It  appears  that  the  larger  the  dP/dt  on  its  rise  to  peak 
pressure,  the  sooner  the  peak  pressure  is  reached.  The  peak  pressure  also 
decreases  but  in  an  insignificant  magnitude  in  these  comparisons.  Figure  4-6 
shows  the  ITP  response  to  varying  dP/dt. 
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Figure  4-5.  Influence  of  abdominal  air  on  ITP  responses 


2.  Isolmpulse  blast  loading  (10.5  psl-ms)  of  square  waves  at  different 
durations  were  also  applied  on  the  surface  of  the  model.  The  higher  the  aagnl- 
tude  of  the  peak  external  blast  and  the  shorter  the  loading  duration  the  fast¬ 
er  the  [d(ITP)]/dt  and  the  higher  that  peak  1TP  occurs.  Also  note  the  peak  1TP 
occurs  at  earlier  tine  (Fig.  4-7).  This  suggests  that  for  the  iapulse  loading 
of  10.5  psl-as  the  higher  peak  external  blast  and  shorter  duration  could  incur 
aore  risk.  Experiaental  results  of  Richmond,  et  al.  (Fig.  4-8)  showed  that  the 
probability  of  lethality  (LD^q)  becomes  much  higher  at  higher  external  blast 
pressure  and  shorter  duration  for  a  given  Iapulse  (dashed  lines).  This  was 
also  demonstrated  In  the  iso-iapulse  study  of  WRAIR  on  sheep. 

4.3.3  Effects  off  Geometric  Rise 

Keeping  the  shape  of  the  model  the  saae  and  retaining  the  saae  material 
properties,  we  vary  the  size  of  the  model  froa  a  "sheep"  to  a  "rabbit"  to 
observe  the  effect  of  animal  size  on  various  physical  parameters  under  the 
same  "short  duration"  blast.  Figure  4-9(a)  shows  different  ITP  responses  due 
to  different  geoaetrlc  sizes.  The  curve  labeled  1.0  corresponds  to  that  of  a 
sheep  considering  the  longitudinal  length  of  the  sheep  torso  to  be  50  ca.  The 
curve  labeled  0.4  corresponds  to  that  of  a  rabbit  with  20  cm  torso  length. 
Comparison  of  the  curves  shows  that  the  smaller  size  animal  is  affected  by  a 
larger  disturbance.  Namely,  the  higher  ITP,  [d(ITP)]/dt  and  higher  oscillation 
frequency.  Figure  4-9(b)  and  (c)  show  comparisons  of  normalized  chest  wall 
displacement  and  front  chest  velocity,  respectively.  Similar  tendencies  are 
seen  on  displacement  and  velocity  histories.  For  a  given  blast  overpressure 
loading  the  small  animal  Is  under  greater  disturbance.  The  ITP  responds  at 
higher  amplitude  and  higher  rate  of  pressure  rise.  The  whole  body  vibrates  at 
a  faster  frequency. 

4.3.4  Effects  of  Varying  Lung  Bulk  Modulus 

By  varying  the  bulk  modulus  of  the  lung  (but  keeping  the  same  Poisson 
ratio,  v  -  0.40)  we  observe  its  effect  on  animal  response  under  the  "short 
duration"’  blase.  Figure  4-10(a)  shows  the  ITP  histories  for  different  cases. 
It  appears  that  for  larger  bulk  modulus  the  ITP  rises  faster  and  a  higher  ITP 
can  be  attained.  The  frequency  response  is  also  higher  under  the  same  external 
blast.  Vo  also  vary  the  bulk  modulus  of  the  abdominal  contents,  rumen, 
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animals  (mice,  rats,  hamsters,  guinea  pigs,  and  rabbits).  All  measurements 
were  made  at  an  ambient  pressure  of  12  psl.  Reproduced  from  Richmond  et  al 
(Ref.  10). 


<•)  ITP 


FE  M  MODEL 


Figure  4-9.  Variation  of  rcaponae  parameter  histories  with  changes  in  geometric 
dimensions  of  the  model*  Sizes  represented  range  from  sheep 
(labeled  1.0)  to  rabbit  (labeled  0.4).  The  geometric  shape  and 
material  properties  are  kept  constant. 
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NORMALIZED  FRONT  CHEST  DISPL  INTRR  THORACIC  PRESSURE  (PSI1 
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FEM  MODEL 


FEM  MODEL 


intestines,  etc*  The  observation  Is  that  It  aakes  a  slight  change  on  ITP  for  a 
wide  range  variation  of  abdoaen  bulk  Modulus*  Figure  4-10(b)  and  (c)  shows  the 
normalised  chest  dlsplaceaent  and  velocity,  respectively*  The  phenomena  is 
similar  to  that  in  ITP  histories*  The  larger  the  bulk  Modulus  of  the  lung  the 
faster  and  higher  the  response  is* 


5.  A  CS0SS—SKCTIOIUL  1HO-ODBMSIOHAL  MODEL 


5.1  nfnooocnoH 

A  three-dimensional  finite  element  model  of  sheep  torso  has  been  con¬ 
structed  as  described  in  the  previous  section.  This  model  16  used  to  predict 
the  intrathoradc  pressure  response  when  a  sheep  is  exposed  to  external  blast 
loadings.  The  prediction  was  compared  with  experimental  measured  esophageal 
pressure  history  as  described  previously.  The  existing  experimental  data  main¬ 
ly  comes  from  overpressure  taken  at  one  or  two  depths  in  the  animal  esophagus. 
Jftnsson  et  al.  [8]  also  measured  lung  overpressures  directly  by  introducing 
transducers  through  the  trachea  into  the  right  and/or  the  left  lung  In  rabbit 
Impact  experiments.  Their  measurements  indicated  that  different  overpressure 
histories  were  experienced  by  different  parts  of  the  lung  in  response  to  the 
Impact  loading  delivered  at  one  side  of  the  body  surface.  Examination  of  the 
lung  trauma  following  incidences  of  blast  waves  or  high  velocity  impact  [8.14] 
reveal  that  gross  hemorrhage  on  the  pleural  surface,  in  particular  the  area 
near  the  Impact  side  end  the  area  near  the  heart.  It  suggests  that  wave 
dynamics  may  play  a  role  in  determining  the  Injury  site  and  severity.  The 
complex  phenomena  of  wave  reflection  at  boundaries  and  wave  focusing  at  geo¬ 
metric  focal  points  may  cause  different  types  of  Injury  at  different  regions. 
To  understand  the  Injury  mechanism  associated  with  the  lung  under  blast  expo¬ 
sure  we  need  to  follow  the  incident  wave  as  it  propagates  across  the  body.  We 
need  to  know  how  different  parts  of  the  body  respond  to  external  disturbance 
and  how  mechanical  stress  or  energy  is  built  up  at  different  regions  of  the 
lung.  With  this  knowledge  wo  can  proceed  to  correlate  the  injury  map  and 
Injury  severity  from  animal  experiments.  An  Injury  criterion  is  expected  to  be 
established  in  terns  of  local  tissue  mechanical  response  parameters. 

To  supply  this  greater  detail,  a  two-dimensional  FEM  model  was  construct¬ 
ed,  This  model  is  complementary  to  the  three-dimensional  model  in  providing 
both  temporal  and  spatial  resolution.  Characteristics  of  wave  propagation 
across  different  parts  of  the  body  can  be  monitored. 
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5.2  DBSCK1PTI0N  CP  THE  MODEL 


5.2.1  Qeo— try 

A  two-dimensional  model  is  still  a  simplification  of  the  actual  three- 
dimensional  body.  To  justify  the  two-dimensional  model  it  has  been  demon¬ 
strated  by  Clemedson  and  Jbnsson  [2,3],  by  using  protective  layers  covering 
either  chest  or  abdomen  (rabbit),  that  lung  injuries  are  caused  primarily  by 
chest  loading.  The  effect  of  mouth  airflow  and  the  Indirect  abdominal  loading 
through  the  diaphragm  were  found  to  be  insignificant  for  lung  injuries.  Based 
on  the  sheep  thorax  cross  section  at  the  level  of  the  seventh  thoracic  verte¬ 
bra  [9],  a  two-dimensional  plane  strain  model  was  constructed  (see  Fig.  5-1). 
The  model  consists  of  four  major  organs  -  skeletal  muscle,  rib,  lung  and 
heart.  The  objectives  are  (1)  to  understand  the  underlying  mechanics  of  wave 
dynamics  in  biological  bodies  during  blast  exposures  and  (2)  to  predict  the 
overpressure  occurrences  at  different  parts  of  the  lung  when  a  sheep  is  under 
blast  wave  incidence. 

Following  the  FEM  formulation  of  Section  3.1,  the  two-dimensional  dis¬ 
cretization  mesh  is  defined  by  133  nodes  and  122  four-node  isoparametric  con¬ 
tinuum  elements.  TWo  degrees  of  freedom  (displacement:)  are  assigned  to  each 
node.  Bilinear  interpolation  functions  are  used  for  element  approximation. 
Implicit  integration  scheme  with  Newmark  approximation  is  employed  for  the 
temporal  discretization.  Analysis  is  performed  for  30  ms  from  the  onset  of 
wave  arrival  at  the  body  surface,  since  experimental  measured  esophageal 
pressures  become  diminishing  after  30  ms  in  response  to  the  loadings  of  inter¬ 
est  in  this  program.  For  most  of  the  analysis  a  time  step  of  0.5  ms  is  chosen. 

5.2.2  Modeling  of  the  Rib  Cage 

The  geometrical  structure  and  the  kinematic  motion  of  the  rib  cage  is 
three-dimensional  by  nature,  so  we  must  find  a  two-dimensional  equivalent 
representation  to  account  for  its  rigidity  and  kinematic  degrees  of  freedom. 
Considering  its  role  in  resisting  as  well  as  delivering  the  external  impluse 
loading  from  body  surface  to  the  internal  organs,  it  is  obvious  that  a  closed- 
ring  of  bony  rib  overexaggerates  its  rigidity.  Yet  a  closed-ring  of  muscle 
also  underestimates  the  global  stiffness  of  the  whole  body  structure.  In  order 
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to  understand  the  Influence  of  rib  cage  representation  on  the  ITP  response  to 
blast  wave  a  parametric  study  with  different  rib  cage  representation  has  been 
carried  out. 

Figure  5-2  (a)  and  (b)  shows  the  closed-ring  and  segmental  type  rib 
representations,  respectively.  For  the  segmental  type  rib  cage  three  kinds  of 
materials  are  used,  namely,  bony  rib,  muscle,  and  cartilage.  The  elastic 
moduli  of  the  cartilage  are  hypothetically  chosen  to  be  ten  times  that  of 
muscle.  Figure  5-3(a)  shows  the  intrathoracic  responses  for  various  rib  cage 
representations.  Curve  1  is  the  response  with  closed-ring  bony  rib.  Curve  2  is 
the  case  with  segmental  rib  composed  of  muscle  and  bone.  Curve  3  is  the  case 
with  closed-ring  representation,  but  with  cartilage  as  material.  Case  4  is  of 
segmental  representation  of  muscle  and  cartilage.  Case  5  is  a  closed-ring 
representation  of  muscle  only.  It  appears  that,  except  for  curve  1,  all  the 
rest  respond  in  similar  frequencies,  but  with  varying  amplitudes.  The  response 
of  case  1  Is  not  reason  ble  due  to  the  overexaggerated  rib  cage  strength. 
Among  the  rest,  case  5  responds  most  similarly  to  experiments  measured  with 
clear  sequential  oscillation  (compression  ♦  tension  +  compression  at  reducing 
amplitude).  It  appears  that  modeling  with  bone  material  deviates  the  ITP 
prediction  from  experimental  measurements.  This  could  imply  that  at  an  early 
transient  period  most  of  the  wave  loading  entered  the  lung  through  muscle 
layers. 

The  corresponding  displacement  histories  of  the  chest  wall  (point  A  of 
Figure  5-2)  is  also  shown  as  Figure  5-3 (b).  The  inward  displacement  of  case  5 
reaches  0.7  cm  at  12  ms.  For  cases  2,  3,  and  4  the  displacements  fall  in  the 
range  of  0.4-0. 6  cm  at  7-10  ms.  For  case  1  the  closed  ring  rib  structure  Is 
equivalent  to  a  very  effective  protective  shell.  Its  rigidity  probably  well 
protects  the  lung  inside.  Therefore,  the  lung  appears  to  be  only  slightly 
disturbed  as  compared  to  other  cases. 

5.2.3  Material  Constants 

Measurements  of  esophageal  pressure  responses  to  blast  wave  loading 
Indicates  that  there  is  strong  damping.  It  first  rises  to  a  peak  pressure 
value  (compression)  then  falls  to  a  negative  phase  (tension)  with  much  smaller 
amplitude  and  Increased  period  time  duration.  Viscoelastic  material 
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description  Is  one  cf  the  damping  mechanisms  causing  such  responses.  The  2-D 
model  consists  of  four  major  organs.  For  each  organ  material,  to  describe  Its 
stress-strain  relationship,  however,  six  constants  are  needed  even  for  the 
simplest  representation.  They  can  be  written  [by  choosing  NG  ■  NK  ■  1  in  Bqs. 
(21)  and  (22)  of  Sec.  3]  as: 


K(t)  -  Kp  +  Kj  exp(-t/A) 

(1) 

G(t)  -  GQ  +  Gj  exp(-t/B)  . 

(2) 

In  Eq.  (1),  K(t)  represents  the  bulk  modulus  at  the  current  time  t,  Kq  repre¬ 
sents  that  when  time  becomes  Infinity,  and  exp(-t/A)  denotes  the  exponen¬ 
tially  relaxing  part  of  the  bulk  modulus.  For  a  given  strain  applied  at 
normalized  time  (t/A)  -  0,  it  is  Kj.  This  terra  approaches  0  as  normalized  time 
(t/A)  becomes  Infinite.  Similar  notations  are  used  for  the  shear  modulus  of 
Eq.  (2).  For  the  four  major  organs  used  in  the  2D  model,  some  of  these  quanti¬ 
ties  have  been  measured,  but  many  have  not.  Before  undertaking  new  experiments 
to  determine  the  missing  quantities,  we  ■will  parametrically  vary  them  to  see 
which  are  Important. 

5.2*4  Senaltlvlty  Studies  of  the  Material  Constanta 

The  major  purpose  of  these  sensitivity  studies  Is  to  determine,  among  the 
total  of  28  material  constants,  which  ones  are  important  and  which  ones  are 
insensitive.  Through  these  studies  the  number  of  material  constants  is  reduced 
to  the  minimum.  Only  the  important  parameters  will  be  retained  and  their 
influences  justified. 

To  make  these  studies  the  case  of  double  peak  blast  loading  with  At  - 
7.6  ms  (time  duration  between  two  blast  peaks)  is  chosen.  It  is  chosen  because 
double  peak  loading  provides  a  sensitive  test  on  the  model.  This  test  allows 
us  t"'  justify  the  model  response  sensitivity  in  both  spatial  and  temporal 
sen  iese  comparisons,  however,  are  relative  in  the  sense  that  the  muscle, 
bone,  and  heart  are  assumed  to  have  the  same  material  constants. 

Viscoelastic  Description  of  Ifasde,  Bone,  and  Heart  -  Figure  5-4  shows 
the  two  model  ITP  response  curves  under  the  given  double  peak  blast  loading. 
The  two  very  close  response  curves  show  that  the  model  is  insensitive  to  the 


Figure  5-4.  I  TP  response  curves  for  esses  with  and  without  viscoelastic 
teras  for  aiscle,  bone,  and  heart. 
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viscoelastic  description  of  the  three  surrounding  organs.  Compared  to  the  lung 
material,  these  three  have  bulk  moduli  close  to  that  of  water.  Whereas  the 
lung  has  a  bulk  modulus  close  to  that  of  air.  These  three  materials  are  Incom¬ 
pressible  when  compared  to  the  lung.  Therefore,  contribution  from  the  visco¬ 
elastic  description  of  these  three  surrounding  organs  In  the  decaying  feature 
of  the  lung  1TP  response  could  be  minimal. 

Shear  Moduli  of  the  Ifascle,  Bone,  and  Heart  -  Magnitude  of  shear  moduli 
of  these  three  organs  are  varied  to  compare  the  effect  on  the  model  1TP 
response.  It  appears  to  be  an  Important  parameter.  Figure  5-5  shows  these 
comparisons.  Among  these  comparisons,  bulk  moduli  Kq  are  kept  at  1.0  x  10 10 
dyne/cm^.  The  case  with  Gq  ■  1.0  x  10*  dyne/ cm ^  is  the  baseline  case.  The  case 
with  Gq  -  1.0  x  10  dyne/cm  shows  little  difference  from  the  baseline  case. 
But  as  we  keep  Increasing  Gq,  the  difference  In  the  model  ITP  response  becomes 
evldlent.  These  cases  show  smaller  amplitude  and  increased  frequency.  This  is 
due  to  the  model  becoming  less  flexible. 

Balk  Moduli  of  the  lfasde.  Bone,  and  Heart  -  Magnitude  of  bulk  moduli  of 
these  three  organs  are  varied  to  understand  their  effects  on  the  model  ITP 
responses.  In  this  study  shear  moduli  of  these  three  organs  are  kept  constant 
at  1.0  x  10*  dyne/cm^.  Bulk  moduli  of  Kq  ■  1.0  x  10*®  dyne/cm^  are  used  as  the 
baseline  case.  There  Is  no  significant  change  In  the  ITP  response  as  we  reduce 

O  O 

the  bulk  moduli  down  to  1.0  x  10  dyne/cm  (see  Fig.  5-6).  Further  reduction 
of  Kq  results  in  slower  ITP  responses  and  larger  amplitudes.  The  extreme  case 
with  Kq  -  1.0  x  10  dyne/cm  represents  the  case  when  bulk  moduli  of  these 
three  organs  are  of  the  same  size  as  that  of  the  lung. 

Bulk  Modulus  of  the  Lung  -  Bulk  modulus  Kq  of  the  lung  is  varied  to  study 
Its  influence  on  the  ITP  response,  while  the  shear  modulus  Gq  is  kept  constant 
at  1.0  x  10*  dyne/cm^.  For  muscle,  bone,  and  heart  the  baseline  values  of 
1.0  x  io*®  and  1.0  x  10*  dyne/cm  are  used  for  the  bulk  and  shear  moduli, 
respectively.  Figure  5-7  demonstrates  these  changes.  It  appears  that  bulk 
modulus  of  the  lung  Kq  is  a  crucial  parameter.  Decrease  by  one  order  of  magni¬ 
tude  results  in  smaller  amplitude  and  slower  starting  time  in  the  ITP  response 
curve,  while  increase  by  one  order  of  magnitude  results  in  very  pronounced 
amplitude  and  earlier  starting  time  in  the  ITP  response  curve.  Increased  bulk 
modulus  means  a  stiffer  lung  with  higher  wave  propagation  speed  (c  -  /K/p  for 
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Figure  5-5.  I  TP  response  curves  at  varying  shear  aodull  for  ausde, 
bone,  and  heart. 
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Figure  5-6*  I  TP  reaponae  curvoa  at  varying  bulk  aoduli  for  wade, 
bone,  and  heart. 
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compression  wave)*  Under  a  given  disturbance,  the  lung  responds  in  a  more 
vigorous  way,  while  the  decreased  bulk  modulus  results  in  an  opposite  v  . 

▼iscoelaatlc  Description  of  the  Lang  -  Experimental  esophageal  pressure 
measurement  shows  strong  damping  in  the  lung  overpressure  response  curve 
during  air  blast  exposure*  Damping  came  from  different  sources*  From  the 
material  point  of  view,  the  viscoelastic  nature  of  the  lung  material  could 
contribute  to  the  damping  of  lung  overpressure  response* 

The  specific  viscoelastic  description  of  the  lung  under  dynamic  blast 
wave  loading  is  unknown  at  this  time.  Assumed  values  are  used  for  the 
relaxation  moduli  to  understand  their  effect  on  the  ITP  response  of  the  lung 
under  dynamic  blast  wave  loading.  Unlike  the  muscle,  rib,  and  heart,  the 
viscoelastic  description  of  the  lung  has  been  shown  to  reduce  the  peak  ITP  and 
the  rarefaction  of  the  response  that  follows  (see  Fig*  5-8).  There  may  be 
other  damping  mechanisms  associated  with  the  intrathoradc  pressure  responses 
to  blast  wave  incidence. 

Effective  Mmdejfttb  and  Ita  Calibration  by  Static  Loadlnit  Experlmenta  - 
The  structure  and  geometry  of  the  lung  and  rib  cage  are  three-dimensional  in 
nature.  Under  blast  loadings  the  kinematics  and  deformation  dynamics  of  the 
rib  cage  and  the  body  are  expected  to  be  three-dimensional  as  well.  To  account 
for  the  structural  rigidity  of  the  rib  cage  in  a  two-dimensional  model  we  have 
lumped  its  rigidity  to  the  material  property  of  the  muscle  layer  and  call  it 
the  effective  muscle/rib.  The  value  of  this  effective  property  can  be  obtained 
from  static  loading  experiments. 

At  the  present  time  a  system  parameter  identification  scheme  is  used. 

Within  the  context  of  the  sensitivity  study,  various  combinations  of  bulk  and 

shear  moduli  have  been  tried  to  determine  the  set  of  parameters  which  makes 

the  best  ITP  prediction  when  compared  with  a  group  of  experimental 

measurements.  In  other  words,  the  mechanical  description  for  the  effective 

muscle/rib  is  obtained  through  a  series  of  numerical  experiments.  Shear 

7  o  q  o 

modulus  of  1*0  x  10  dyne/cm  and  bulk  modulus  of  1*0  *  10  dyne/cm  have  been 

found  to  be  the  pair  from  this  comparison  study.  Their  values,  in  conjunction 

with  other  elastic  properties,  are  listed  in  Table  5-1. 
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Figure  5—8.  I  TP  response  curves  for  esses  with  end  without  viscoelastic 
terms  for  the  lung* 
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Table  5-1*  Material  Properties  Used  In  Cross-Section  Model 


Organ 

°0 

(dyne/cm2) 

*0 

(dyne/cm2) 

P 

(gm/cm3) 

Lung 

1.0  x  10A 

1.4  x  io6 

0.2 

Muscle 

1.0  x  107 

1.0  x  io9 

1.0 

Rib 

1.0  x  107 

Os 

o 

w 

X 

o 

• 

2.0 

Heart 

1.0  x  106 

1.0  x  1010 

1.0 

Static  loading 

experiments 

can  be  done  to 

confirm  the  effective 

auscle/rlb  description.  At  this  tine,  preliminary  work  is  done  on  the  model  to 
justify  this  type  of  experiment.  Static  loading  Is  applied  on  the  model  chest 
wall  with  the  back  constrained  and  the  chest  wall  deflection  is  predicted. 

Figure  5-9  shows  the  body  configurations  of  the  unloaded  and  loaded 
states  at  various  degrees  of  static  pressure.  The  deflection  of  point  A  under 
different  static  loadngs  is  plotted  in  Figure  5-10(a).  The  normalised  deflec¬ 
tion  of  point  A  (with  respect  to  the  cross-section  dimension  of  AB)  Is  plotted 
In  Figure  5-10(b).  The  result  shows  that  for  lew  pressure  static  loading 
(below  4  psl)  the  deflection  Is  In  a  reasonable  range. 

5.3  OOMPAKISOH  WITH  KXPKKIMKHTAL  DATA 
5.3.1  Blast  loading  Description 

When  an  Incident  blast  wave  front  (compression)  reaches  the  body  surface 
a  sequence  of  subtle  Interactions  can  follow.  Depending  on  the  relative  impe¬ 
dance,  a  portion  of  the  wave  energy  transmits  into  the  body  and  the  other 
portion  reflects  from  the  interface.  The  reflected  component  is  added  to  the 
newly  arriving  blast  wave  to  result  in  increased  loading  strength.  At  the  same 
time  different  types  of  waves  (shear,  surface)  can  be  set  into  motion.  In 
accordance  with  the  dynamic  equilibrium  among  external  forces,  system  iner¬ 
tial,  elastic  and  damping  forces  the  body  deforms  and  the  body  surface  moves. 
At  a  given  instant,  the  blast  loading  on  the  body  surface,  in  fact,  depends  on 
the  entire  past  blast-body  interaction  history. 


(e)  4.6  pal  (d)  6.9  pal 

Figure  5-9.  Deformed  configuration  of  the  two-dimensional  FEM  model  under 
prescribed  static  loadings* 
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The  blast  loading  description  Is  best  described  in  terms  of  an  Impulse 
ratio.  The  Impulse  ratio  is  defined  to  be  the  factor  by  which  impulse  (time 
integration  of  pressure  history)  is  increased  at  the  body  surface  over  that  of 
the  free  field  value.  Based  on  the  experimental  data  reported  in  Reference  20, 
it  is  assumed  that  the  impulse  ratio  at  the  body  surface  can  be  expressed  as  a 
function  of  external  peak  blast  pressure.  Figure  5-1 1(a)  illustrates  the 
derived  values  from  experimental  data  and  the  hypothetically  assumed  linear 
relationships  between  impulse  ratios  and  external  peak  blasts.  Data  from  the 
experiments  are  marked  with  X.  The  hypothetical  linear  relationship  is 
obtained  by  extrapolating  the  two  experimental  values.  Division  of  zones  A,  B 
and  C  are  illustrated  in  Figure  5-ll(b). 

Zone  A  is  the  portion  of  the  body  surface  directly  facing  to  the  incident 
wave.  Between  0  and  45  psl  blast  peak  pressure  the  impulse  ratio  varies 
linearly  from  1  up  to  2.8.  The  increased  impulse  ratio  accounts  for  the  effect 
of  wave  reflection  from  body  surface  and  the  resulting  higher  loading*  Zone  C 
is  at  the  back  side  of  the  body  away  from  the  direction  of  wave  arrival.  When 
the  free  field  blast  pressure  is  3  psi  the  impulse  ratio  derived  from  experi¬ 
mental  data  is  found  to  be  0.75.  In  Zone  C  it  is  assumed  that  the  ratio  is 
0.75  for  peak  pressures  less  than  10  psi  and  1  for  blast  pressure  higher  than 
10  psi.  Zone  B  is  the  body  surface  normal  to  the  direction  of  wave  arrival. 
The  impulse  ratio  in  Zone  B  is  assumed  to  be  the  average  values  between  Zone  A 
and  Zone  C.  This  body  surface  blast  loading  is  used  in  all  the  studies  in  this 
report . 


5.3.2  Spatial  Variation  of  Overpressure  Histories 

Figure  5-12(a)  is  reproduced  from  experimental  results  performed  at  the 
Lovelace  ITRI.  A  sheep  is  exposed  right  side  on  to  blast  waves  generated  with 
32  lb.  of  TNT  charge  at  44.4  ft  range  and  8.0  ft  height  of  burst.  The  free 
field  blast  pressure  and  the  intraesophageal  pressure  taken  at  the  seventh 
intercostal  are  shown  as  chain-dotted  and  chain-dashed  lines,  respectively. 

Figure  5-12(b)  shows  the  resulting  lung  overpressure  responses  as  pre¬ 
dicted  by  the  two-dimensional  sheep  thorax  model.  Each  of  the  seven  curves 
denotes  ovevpressure  history  occurred  at  a  different  location  in  the  lung  [see 
Fig.  5-ll(b)jo  It  is  clear  that  the  incident  wave  arrives  at  different 
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(a)  Effect  of  experimental  derived  and  hypothetical  Impulse  ratios 
at  various  peak  blast  pressures  on  Zones  A,  B  and  C 
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(b)  Division  of  body  surface  Into  Zones  A,  B  and  C 
with  respect  to  the  explosion  center 


Figure  5-1 1.  Experimental  derived  and  hypothetical  impulse  ratios 
at  various  peak  blast  pressures. 
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Figure  5-12.  Experimental  free  field  blast  pressure,  esophageal  pressure  and  overpressures 
predicted  by  the  two-dimensional  FEM  model. 


stations  at  different  times*  The  overpressure  experienced  by  different  parts 
of  the  lung  is  also  different  in  its  temporal  variation  on  account  of  the 
irregular  geometry,  complex  structure,  as  well  as  the  dynamical  behavior*  It 
is  seen  that  station  1  is  the  first  to  "feel"  the  disturbance  followed  by 
station  2.  The  station  1  is  able  to  experience  a  higher  overpressure  and  a 
wider  time  span  due  to  the  fact  that  higher  Impulse  is  delivered  to  its  neigh¬ 
borhood*  Station  3  is  the  third  to  "feel”  the  incident  wave*  Compared  with 
station  1,  station  3  experiences  a  wider  time  span,  but  lower  peak  pressure, 
although  both  are  near  the  right  body  surface  under  Increased  loading.  This  is 
because  station  1  is  located  in  a  narrow  corner  between  skeletal  muscle  and 
the  heart.  When  the  reflected  component  from  the  heart  is  superimposed  to  the 
Incident  wave,  a  higher  overpressure  can  occur  at  this  region.  Most  of  the 
wave  passing  by  station  3  Is,  however,  able  to  penetrate  into  the  lung  easier* 

Station  4  is  the  fourth,  among  the  seven,  to  know  the  pressure  wave  inci¬ 
dence.  It  rises  to  1  psi  at  4  ms  yet  experiences  a  second  peak  of  4  pel  at  10 
ms.  This  second  peak  is  caused  primarily  by  the  arrival  of  the  pressure  wave 
which  already  passed  stations  3  and  7.  Station  5  is  the  fifth  to  "know"  due  to 
the  fact  that  it  is  deeper  from  the  body  surface  and  also  has  a  lower  pressure 
applied  at  that  surface.  Station  6,  which  is  located  behind  the  heart,  is  the 
next  to  respond*  Station  7,  at  a  location  approximating  the  center  of  the  lung 
where  the  esophagus  sits  (we  have  neglected  the  mediastinum,  blood  vessel,  and 
soft  tissue  there),  is  the  last  one  to  know  the  pressure  arrival  yet  reaches  a 
maximum  value  of  6  psi  at  7  ms. 

Figure  5-12(a)  also  shows  the  intraesophageal  pressure  history  taken  at 
the  seventh  Intercostal  from  sheep  experiments.  Compare  with  this  experimental 
data,  all  prediction  curves  at  seven  stations  show  higher  frequencies,  but 
similar  average  response  with  the  removal  of  high  frequency  components.  To 
make  the  comparison  we  divide  the  first  30  ms  into  three  time  periods  on  both 
Figure  5-12(a)  and  (b): 

1.  Time  period  A:  from  0-12  ms 

2.  Time  period  B:  from  12-25  ms 

3.  Time  period  C:  from  25-30  ms 
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It  Is  seen  that  during  the  time  period  A  all  the  seven  stations  experience 
overpressure  due  to  the  arrival  and  reflection  of  pressure  waves  in  different 
ways  at  different  times.  During  the  time  period  B,  however,  almost  all  of  the 
seven  stations  are  under  different  degrees  of  dilation*  The  physical 
Interpretation  is  that  the  whole  lung  is  progressively  under  recoil  tension 
after  the  compression  In  the  period  A.  The  decreased  amplitudes  and  retarded 
oscillation  frequencies  are  attributed  to  the  damping  in  the  lung.  After  25  ms 
all  the  response  curves  in  the  time  period  C  show  their  trend  to  bounce  back 
to  the  compression  phase  with  less  available  energy.  Compared  with  esophageal 
pressure  history  measured  from  the  experiments,  prediction  curves  of  Fig¬ 
ure  5-12 (b)  show  that  there  is  spatial  differences  of  lung  overpressure 
responses  to  external  blast  incidence.  Locally,  different  response  history 
will  result  in  different  damage  history.  Therefore,  the  interpretation  of 
esophageal  pressure  should  be  carefully  justified  as  it  is  commonly  used  as  an 
indication  of  the  whole  lung  uniform  response  in  experimental  works  as  well  as 
lumped-parameter  model  prediction. 

5.3.3  Single  Beak  Exposure 

Data  of  external  blast  loading  taken  from  WRAIR/ LOVELACE  ITRI  experiments 
are  applied  on  the  two-dimensional  FEM  model  to  predict  lung  overpressure 
responses.  Prediction  of  lung  overpressure  histories  at  the  center  of  the 
lung,  approximately  where  the  esophageal  pressure  is  taken,  are  then  compared 
with  experimental  measurements. 

Figures  5-13  through  5-16  are  comparisons  of  the  four  different  loading 
cases  from  TNT  free  field  experiments  performed  during  September  and  October 
of  1981.  In  each  case  the  model  prediction  is  overlayed  on  the  experimental 
esophageal  pressure  recording  for  easy  comparison.  The  comparisons  show  that 
the  two-dimensional  model  prediction  and  the  experimental  measurements  are  of 
good  qualitative  agreement.  Similar  to  the  case  as  demonstrated  in  Section 
5.3.2  and  in  Figure  5-12,  the  overpressure  responses  are  predicted  to  be 
different  at  different  parts  of  the  lung.  Without  the  availability  of 
experimental  measurements  the  multiple  prediction  tracings  are  not  shown  on 
Figures  3-13  through  5-16. 
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Model  prediction 
Esophageal  overpressure 


Figure  5-13.  (Case  1)  8  lb  TNT  free  field.  Experimental  free  field  blast  pressure,  esophageal 

pressure  and  overpressure  predicted  by  the  two-dimensional  FQI  model. 
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5*3*4  Doable  Pwk  Exposure 


The  double  peek  study  was  performed  by  the  LOVELACE  1TRI  in  October  1981 
and  May  1982*  The  purpose  of  the  study  was  to  investigate  resonance  effects 
when  the  tine  duration  between  two  blast  loading  approaches  that  corresponds 
to  the  natural  frequency  of  the  structure*  Experimental  blast  loading  data 
with  different  tine  duration  between  two  peaks  are  applied  on  the  two- 
dinenslonal  FEM  model  body  surface  to  observe  the  lung  overpressure  responses. 
Figures  5-17*  5-18  and  5-19  are  comparisons  of  the  three  cases  with  duration 
between  two  peak  blasts  At  -  9.7*  7*6  and  3.6  ms*  respectively.  Similar  to  the 
previous  discussion  the  overpressure  histories  at  seven  statons  are  plotted  to 
compare  with  the  experimental  measured  esophageal  pressure. 

In  the  third  case  (Fig.  5-19)  with  At  *  3.6  ms,  the  model  prediction  does 
not  clearly  distinguish  the  two  ITP  peaks.  Further  refinement  of  the  finite 
element  mesh  is  needed.  For  the  first  and  second  cases  with  At  -  9.7  and  7.6 
ms  (Figs.  5-17  and  5-18),  model  predictions  compare  favorably  with 
experimental  measurements.  The  model  predictions  do  not  seem  to  have  enough 
rarefaction  as  experimental  measurements.  Typically*  they  have  larger  tension 
phase  and  continuous  ripples.  This  is  due  to  Insufficient  damping.  This  is  to 
be  improved  in  the  next  phase. 

More  comparisons  of  model  prediction  and  experimental  measurements  are 
reported  in  Section  5.4.2  under  "Isoimpulse  with  Different  Peak  Blast 
Pressure." 

5.4  SEBSITIVITT  STUDIES  OF  THE  MODEL 

5*4*1  Effect  of  Loading  Firing  Sequence  (temporal)  and  Doubling  Effect 
Djatrjlwtloo  (spatial) 

The  Interaction  between  the  incident  blast  wave  and  the  deformable  body 
is  a  rather  complicated  coupling.  The  Incident  blast  wave  front  reaches  dif¬ 
ferent  parts  of  the  body  surface  at  different  times  on  account  of  its  propaga¬ 
tion  in  the  air  as  well  as  the  Interaction  with  the  deformable  body.  The 
reflection  of  the  blast  wave  at  the  boundary  also  results  in  increased  effec¬ 
tive  loading  pressure.  Experimental  measurements  and  numerical  simulations 
have  been  done  with  a  solid  torso  model  to  study  its  interaction  with  the 
Incident  blast  wave  [20]. 
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Figure  5-19.  (Case  21)  Double  peak  study.  At  -  3.6  ms.  Experimental  free  field  blast  pressure 
esophageal  pressure  and  overpressure  predicted  by  the  two-dimensional  FEH  model. 


In  order  to  understand  the  effect  cf,  (1)  increased  loading  factor  spa¬ 
tial  distribution,  and  (2)  temporal  variation  in  the  time  of  wave  arrival  on 
the  animal  response,  the  model  is  divided  into  five  different  zones,  as  illus¬ 
trated  in  Figure  5-20(a).  Directly  facing  the  direction  of  wave  arrival  zone  A 
is  the  first  part  of  the  body  surface  under  blast  wave  incidence  and  is  the 
part  under  the  most  extensive  increased  loading*  Follow  the  order  of  B,  C,  and 
D  different  zone  starts  exposing  to  the  incident  wave  at  a  sequentially 
delayed  time  and  reducing  loading  factor  of  1*8,  1*6  and  1.4,  respectively. 
Body  surface  at  zone  E  is  the  last  to  experience  the  wave  and  has  no  appreci¬ 
able  Increased  loading.  Zone  E  is  assumed  to  lag  zone  A  by  2  ms.  Histories  o' 
increased  loading  at  different  zones  are  depicted  in  Figure  5-20(b).  Figure 
5-20(c)  shows  the  lung  overpressure  at  seven  stations  under  the  prescribed 
loading  distribution  and  timing  sequences.  Compared  with  Figure  5-20(d),  in 
which  the  model  is  under  the  same  spatial  loading  distribution  but  applied  si¬ 
multaneously,  it  is  seen  that  only  station  2  and  station  4,  which  is  located 
in  the  back  side  of  the  lung,  start  to  respond  to  the  incident  wave  at  a 
delayed  time.  No  significant  difference  is  observed  in  the  responses  at  the 
rest  of  the  transducers. 

5*4.2  Iso imp  alee  Stmdlee 

To  further  examine  the  responsive  behavior  of  t.«e  two-dimensional  model 
Isoimpulse  blast  loading  is  applied  to  the  model  in  the  following  three  ways: 

1.  Isoimpulse  with  different  dP/dt  (triangular  waves  of  the  same  dura¬ 
tion,  sample  peak) 

2.  Isoimpulse  with  different  peak  blast  pressure  (isosceles  triangular 
waves) 

3.  Isolmpulse  with  different  loading  frequencies  (half  sine  waves) 

Isoimpulse  with  Different  dP/dt  (triangular  waves  of  the  same  duration, 
sample  peak)  -  A  constant,  impulse  is  applied  on  the  body  surface  at  different 
dP/dt  (triangular  wave  at  fixed  peak  blest  pressure  and  loading  duration).  The 
ITP  responses  at  station  7  are  recorded  for  comparison.  Figures  5-21  ana  5-22 
are  results  for  Impulse  ■  14  and  23  psl-ms,  respectively.  It  appears  that  the 
steeper  the  dP/dt,  the  faster  the  ITP  prediction  rises  up  and  resulting  in  an 


(a)  Division  of  five  different  (b) 
cones  (labeled  as  A,  B,  C, 

D,E). 


Histories  of  Increased  loading  factors 
at  cone  A,  B,  C,  Dt  E,  respectively* 
(Note.  Follow  the  order  of  A,B,C,D,E* 
There  lo  a  reducing  loading  factor  as 
well  as  a  delayed  starting  tiae.) 


Figure  5-20*  Effects  of  varying  firing  sequence  and  loading  factor  histories 
on  the  lung  overpressure  responses* 
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(c)  Resulting  lung  overpressure 
histories  with  varying 
firing  sequence. 
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(d)  Response  curves  with  simultaneous 
loading 


Figure  5-20.  (Cont'd) 
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Figure  5-21.  Isoinpulse  (14  pai-as)  study  with  different  dP/dt 
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earlier  peak  overpressure  cine.  The  magnitudes  of  the  peak  ITPs,  however*  are 
approximately  the  same  In  these  comparisons* 

Isolmpulse  with  Different  Peak  Hast  Pressure  (Isosceles  triangular 
waves)  -  For  short  duration  air  blast  waves*  the  Impulse  has  been  considered 
as  an  Important  parameter  In  causing  body  injury  or  even  death*  For  a  given 
Impulse*  however*  the  resulting  injury  also  depends  on  the  peak  pressure*  A 
fixed  amount  of  impulse  (22*1  psl-ms)  was  used  in  the  WRAIR  sheep  exposure 
experiments.  This  lsolmpulse  is  applied  on  the  animal  with  different  combina¬ 
tions  of  peak  blast  pressure  and  time  course.  The  measured  maximum  I  TP  (the 
esophageal  pressure)  and  the  corresponding  calculated  maximum  d(ITP)/dT  (the 
I  TP  increasing  rate)  are  plotted  against  the  peak  blast  pressure  on  Fig¬ 
ure  5-23(a)  and  (b)*  respectively.  The  resulting  gross  lung  injury  was  divided 
into  three  different  groups  separated  by  two  marks  indicating  borders  of  "no 
lung  injury"  and  "severe  lung  injury". 

In  order  to  understand  tha  effect  of  peak  blast  pressure  and  the  result¬ 
ing  I TP  responses  at  constant  impulse*  an  impulse  of  22.1  psi-ma  is  applied  on 
the  model  at  different  peak  blast  pressure  (isosceles  triangular  waves).  Ten 
different  loadings  whose  peak  blast  pressure  varying  from  2.46  up  to  44.2  psl 
were  applied  on  the  model  and  the  transient  interaction  are  followed.  To 
compare  with  the  WRAIR  experimental  results  under  the  sane  type  of  loading* 
the  maximum  ITP  and  the  maximum  d(ITP)/dT  from  model  prediction  are  plotted  as 
dotted  lines  on  Figure  5-23(a)  and  (b)  of  WRAIR  previous  results  for  compari¬ 
son*  The  model  prediction  at  the  center  of  the  lung  (the  subplots  at  the  upper 
right)  is  taken  for  this  comparison.  Since  it  is  at  a  comparable  location 
where  the  esophageal  pressure  was  measured. 

Detailed  ITP  responses  of  six  representative  cases  are  shown  in  Fig¬ 
ure  5-24(a)  through  (f).  For  each  case  the  upper  left  subplot  shows  the 
external  blast  pressure  history.  The  rest  of  the  nine  subplots  show  the  ITP 
responses  at  various  locations  in  the  lung  as  indicated  by  arrows. 

A  study  of  the  responses  at  various  loading  cases  suggeots  that: 

a.  At  higher  blast  pressure  and  shorter  duration*  the  wave  feature 
becomes  important.  Wave  focusing  at  the  geometrical  focal  point  and 
wave  reflection  from  a  neighboring  "firm  organ"  can  lead  to  local 
nigh  overpressure  response. 
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and  luaped-paraaeter  predictions 
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(b)  Peak  pressure  -  29.5  psl 
Figure  5-24.  (Cont'd). 
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(e)  Peak  pressure  -  9.82  psi 
Figure  5-24.  (Cont'd). 
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b.  Probe  1  is  located  In  the  wave  lapact  aide  at  a  corner  between 
Incompressible  material,  where  wave  reflection  1c  likely  to  occur 
aore  frequently  In  a  given  time.  Repetitive  additive  effect  at  such  a 
narrow  cone  can  result  In  very  pronounced  local  loading*  Probe  2  Is 
also  located  at  a  corner  between  Incompressible  material  but  at  the 
back  side*  Similar  wave  reflection  results  In  high  local  pressure* 

c.  Probe  7  Is  located  approximately  at  a  geometric  focal  point.  The 
results  show  the  wave  focusing  at  this  point* 

d*  At  low  peak  blast  pressure  and  long  duration  a  uniform  response  Is 
expected  In  the  whole  lung  [see  Pig.  5-24(f)]. 

Iaolmmulse  with  Mfferemt  Loading  Fraqwamciem  (half  sine  waves)  - 
Isolmpulse  blast  loadings  of  half  sine  wave  with  varying  frequencies  and  peak 
blast  pressure  are  also  applied  on  the  model  to  see  the  variation  In  I  TP 
predictions*  Three  cases  with  frequencies  equal  to  111,  56,  and  28  Hz  (the 
corresponding  angular  frequencies  u  are  700,  350,  and  175  radlan/sec,  respec¬ 
tively)  were  carried  out,  see  Figure  5-25*  In  the  first  case  with  111  Ha,  the 
transient  response  In  different  parts  of  the  lung  Is  clearly  seen*  However, 
for  the  case  of  28  Hz,  even  with  the  same  amount  of  Impulse  loading,  the  whole 
lung  responds  In  a  nearly  uniform  way.  This  demonstrates  that  when  under  high 
frequency  blast  wave  loading  different  parts  of  the  lung  experience  different 
loading  history  on  account  of  the  wave  propagation  effects*  Por  the  low  fre¬ 
quency  loading,  the  external  forcing  function  has  enough  time  to  become  domi¬ 
nant  to  drive  the  whole  lung  to  respond  In  a  uniform  way*  If  the  local  IT? 
history  la  an  Indication  of  the  mechanical  stress  the  local  lung  tissue 
experiences,  the  former  case  Is  likely  to  create  more  Injury  In  terms  of  area 
(or  volume)  and  degree  of  severity*  And  damages  in  the  latter  case  are  more 
duo  to  gross  body  deformation.  Jbnsson  et  al*  [8]  indicate  that  severe  hemor¬ 
rhages  In  the  lungs  were  produced  at  impact  velocity  exceeding  10  m/s.  Por 
velocities  below  5  a/s  acconpanying  large  body  deformation  lethal  injuries 
were  produced  without  sever*  hemorrhages.  Por  a  given  Impulse  short  duration 
high  peak  blast  Is  likely  to  create  higher  chest  wall  velocity  and  result  in 
damages  of  the  former  category*  An  Injury  map  In  terms  of  area,  penetration 
depth  from  pleural  surface,  and  degree  of  severity  with  respect  to  geometric 
descriptions  [curvature  of  the  local  pleural  surface,  bordering  organs  (ribs, 
heart,  etc*)]  and  anatomical  descriptions  (hilar  support,  large  vessel, 
airway,  etc*)  will  assist  In  the  assessment  of  damage  mechanism. 


SHEEP  HOTEL  SHEEP  HOTEL 


Figure  5-25.  Isolmpulse  study  with  varying  loading  frequencies  and  anplitudes 
Case  1:  (a)  &  (b).  Case  2:  (c)  ft  (d),  Case  3:  (e)  &  (f). 


5*4.3  Effect  of  Fat 


In  order  to  understand  the  effect  of  fat  on  the  ITP  responses  when  a 
sheep  Is  under  blast  exposure,  an  additional  layer  representing  fat  is  added 
to  the  body  surface  cf  the  sheep  model*  The  material  is  represented  by  an 
incompressible  gel  type  material  with  density  and  bulk  modulus  the  same  as 
water*  The  thickness  of  the  fat  layer  is  increased  up  to  6  cm  and  the  corres¬ 
ponding  IT?  histories  at  the  center  of  the  lung  region  are  shown  on  Figure 
5-26. 

Due  to  the  increased  total  weight,  the  added  fat  layer  appears  to  reduce 
the  peak  ITP  as  well  as  to  slow  the  ITP  response.  However,  the  degree  of 
change  (8Z  decrease  in  peak  ITP  for  6  cm  increase  in  fat  layer)  does  not  seem 
to  be  very  significant  within  the  range  of  increase. 
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Figure  5-26*  The  model  ITP  responses  (at  the  center  of  lung  region)  at 
varying  fat  thicknesses* 
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A  simplified  16-element  three-dimensional  model  is  constructed  in  Section 
4.  This  nod  el  is  shown  to  reproduce  the  major  body  response  of  a  sheep  under 
blast  wave  loading.  To  study  the  wave  propagation  characteristics  inside  the 
thorax  in  greater  detail,  a  two-dimensional  model  is  presented  in  Section  5. 
This  aodel  shows  that  different  regions  of  the  lung  can  experience  different 
overpressure  histories  during  exposure.  It  suggests  that  the  interpretation  of 
measured  esophageal  pressure  should  be  carefully  justified  as  progressive  wave 
propagation  could  result  in  different  response  histories  at  different  parts  of 
the  lung  and  thereby  different  damage  mechanisms. 

The  model  is  validated  with  experimental  results  from  various  loading 
cases,  e.g.,  4  single-peak,  3  double-peak,  and  10  cases  from  lsolmpulse 
studies.  Reasonably  good  agreement  is  obtained  among  comparisons  made  between 
the  model  I  TP  predictions  and  the  experimental  esophageal  pressure  measure¬ 
ment. 

Material  parametric  study  is  carried  out  to  identify  the  material  parame¬ 
ters  which  are  important  in  the  body's  ITP  response  to  external  blast  wave. 
These  Include  the  density,  the  shear  and  bulk  moduli  of  the  effective 
muscle/rlb,  and  the  bulk  modulus  and  the  density  of  the  lung. 

More  refinement  and  Improvement  on  the  current  model  is  expected  in  the 
next  phase.  Refined  spatial  resolution  will  certainly  reduce  the  numerical 
error  in  the  solutions.  At  high  bl^  ,*t  loading  pressure  (ppeafc  >  25  psi)  the 
ITP  response  does  not  have  as  high  peak  pressure  as  the  experimental  measure¬ 
ment*  This  is  to  be  Improved  by  using  the  ideal  gas-type  nonlinear  stress- 
strain  law  of  the  lung.  The  insufficient  rarefaction  of  the  ITP  responses, 
especially  that  at  the  tension  side,  can  be  improved  by  proper  treatment  of 
damping  mechanisms.  These  Improvements  are  tasked  in  the  next  phase. 
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OIL.  BAMO/UC11 


UTLC!  ).R£AD8S/ASClI2?lfkl 


Identification  division. 

PROfA^-TU.  PtALTT._ 

ENVIRONMENT  DIVISION. 

CONFIGURATION  SECTION. 

SOURCC-CQHPUTiP.  l^IVAr-UHQ. .  . . . 

06JECT-C0KPUTEP.  'jw; VAC- 11^2. 

Ii.PuT-rdTPUT  SrLT*CN. 

_ riLC-Cr\TROL.  _ _  _  _ _ 

SELECT  InFILE  ASSIGN  TO  INTERCHANGE  Y • 

SELECT  OUT  ASSIGN  TO  PRINTER  2. 

DATA  DIVISION. _ _ _ _  _ 

FILE  SECTION. 

FO  INFILE  L*i.EL  RECOkO  IS  OMITTED  RECORDING  MODE  F 
.  tf LOC|L„C on  t  a  t NS.  . .  ...  ,;a  .  RECORDS.. . 

?1  I NR EC  ?iC  X  A  ® ^  . 

ro  OUT  LAGEl  nECOR'P  IS  OMITTED. 

?X  cwtr^c..  Pic.xt.et>.* _ _ _  _  _  . 

i’0ftKZN6  *ST0RA  wE  SECTION. 

?1  C Of UTI T  »:c  9SR9R  VALUE  J. 

PROCEDURE  riVlSION. 

s’tartit. 

OPEN  INPUT  INFILE  »ITH  NO  REWIND. 

.  OPEN  OUTPUT  JOUT .  . _  .  _  _ 

LOOP. 

"EA3  I*  FTuE  AT  END  00  T0  THEND . 

ADD  i  TO  COUNT JT . 

*6Vt  IN^EC  TO  CUT^EC. 

VRITE  OUTAEC. 


GO  TO  LOOP. 
ThE.V). 


close  inf;le 

*TTH 

'‘LOSE  OUT. 

display  * 

COwNT 

STOP  Pu^. 

pc  E  *  I N  G  . 
COUNTIT 


UPON  PRINTER. 
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»  o  i  n  o.i  i  n  r>  o  i  /» 


LiK.Nonmu 


LT'h!,  >.*UVFEA1  I  U) 

P*OGeA,J  **  0  7  F  E  A 

this  PPGORAm  HEADS  OUTPUT  (TAPE  lb)  GENERATED  PU  NCkUTlLlT*  of 
MOVIE. *» VU  .  DATA  ON  TAPE  18  IS  IK  MOVIE  FOK* AT .  PROGRAM  KOVFEA 
WILL  REARRANGE  THE  J A T A  IN  FEAP  FORMAT  AND  W&ITE  IT  TO  T APF  19. 
ICPD  A®«A V  IV  MODaL  CLOCK  IS  TO  9E  SUPPLIED  MC»E . 

“AT  A*»HAy  IN  ^urME^',■  CLOCK  IS  TO  PE  SUPPLIED  MEPt. 

THIS  P’OGRAN  P»>LPapl$  The  FOLLOWING  3  rLQGKS  i 
FEm^-P 

ELEMt NT 

PV’A»VETEP"  NUMMA*P=»  T  f.lJHVPPsis" . »  MJHvLP=lV 

DI^EVSTCN  NPU’.MUMHAPI,  X  ( I  ,  f.UMNPP  )  ,  IP ( 6 , NUHELP ) , 

.  .  _ I  CO D  ( UUMKP.P  )_*___  H. *  T  j K  UHE  l p » _ 

-  !*EAO  INPUT  FPOM  NEWUTILITY  - 

9E*G  US  Vi  S^rVjMKAT.Nim'NPtNUH^LT 
*»UWMATrIA6l(NU?<.sAT) 

PLACf  1®  ,1  5*)  (  I'JfLil.  J)  ,1  =  1.2) ,  J  =  !  ,NJNMAT) 

pEAoa®,  laEiuxti.Ji  ,:=i,i) ,  j=i,nu«np) 

RE«DC19,15D  J IIP (I, J), 1=1,®) ,J=1,NUMELT) 

53  FOsK A T 1 16  IS ) 

ST  ro’-KATCbU?.?) 

-  vs  nr  :uittali_7ation  capos  - 

NLA's  *8 
’ LET'HrG 
i ► rc =1 
lJV*N  =  l?- 

WRITE! 1 9, *) *  fLAK’-'R  -  TEST  FUN  OF  LIKK.MGVFEA  -  ' 

WRITEI13,*)*  3 x  Y 2* 

WRITE! 1«,*) *  "  3  U  V  k* 

WRITE!  19,  i31DNuMNP,  VUMELT,I.UKMAT,NEN,NfcXTPA,  IPEC.^CAN 
UID  F QRhAT  !  71  S ) 
r 

C  -  WRIT**  NOE-AL  CAROS  - 

TCCO!  l)S3!r;3 
TCOO !  2  )  =  -UCO** 

ICOO!  TjroriCOO 
ICO D!  #)=;?i*D3 
ICCDf'f  )  =  y»-P'uE' 

TCOD!  £)=C:CCD3 
ICOO!  ®)  =  D!C“jr 
ICODflTlsSl&rC?  ' 

ICOO ! 13 )  =  j! 1CD3 
ICOGI 19  »  =  30i20S 

‘  icoc  1 1*  )':&i;ri?  ~  . . .  . . . 

icod  i  i"M  =  :u  cr  jO 

:c"Gi i®)=j*icj* 

ICOO*  i"’  )  =  .lnlDD3 
ICCD12E )=uOlCDr 
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i  coo  1 2M=o<ncsn 
icooii^is ' 
icoot3»):uri  ra? 

TCOO( 371=011000 

TCCtM  3* >  =  01lOCO . . . 

TC0DI3cl=01DC03 
ic^ou'is  jur«ir 
iceo««*i  1=  :i23o? 

TCO0l4?)  =  01G"0f!  _ 

lC0D»'«»!f=3!iC0O 
ICCOtH«Or  jt  orcc 
ItC0<4T)= 

WHITE  I !>♦*)•  ‘.'wPaL’ 

00  i‘0 

K*3  WHITE  ns*  irsOJ.IwOuU),  CXII,JI»T=1,7> 
WRIT?  fj ?» ) 

1320  FO^MATtlE  *IlSf?F)w*3) 

1021  F0OM«TnH_l _ _ 

r  . . 

r.  -  WRITE  ELEMENT  EA&DS  - 

"AT  (  1U1  _ _ 

"AT(  ii=l 
"AM  31=1 

"ATJ  U!=1 
A T  (  El  l? 

"AT (  o!=? 

"AT(  71  =  7 
>AT'i  S)=T 
’•AT  (  91=4 

t*AT(  lCI=“  __  _ _ 

MAT'aii=r 

MATI12»S* 

"AT(  1 31  =  *• _ _ _ _ 

»*ATn4>s€  * . . 

M*Ttlil=6 
•*ATt  1*1=* 

vHTT?nst#i •  '" "  rLri*rNT  * 

00  :“Ci  jSlfMIM?LT 
VHITEU9,  i*3pl  J,^AT|  J|  , 

733  WRTTtnS,  10511  U>t  I,  J)  ,i  =  i,6) 

VRTTEC19, 1021 1 
1330  F03MAT<4IS»iri3> 

!  I  J1 F  jCM/TT'f  Z?  14  ) 
r 

STOP 

END' 
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LXMK.nAWV 


LI»>ti*.rt«KOvCin 

C  PROGRAM  FEtHOV 

C  TM£  PROGRAM  fTA^OV  _PE.AUS  'HJTPUT  ION  TAP?  ?r>t  E,6%  DISPLACEMENT 

C  AMO  STRESS  COMPONENTS  AT  EACH  TIME  STEP,'  FROM  FEAP  AFTER  FEM 

C  TPANSIEVT  ANALYSIS  AND  WRITES  VECTOR  ANP/OR  scalar  functions 

c  im  mov if  format  to  tape  n. 

c 

niMEMSTOM  XY?C  ■»,«•->»  U(  T,4S>,X1G 17,7,16) 
«fE8Di:*,UruH.TSlNr>lMtMLFtHU“NP,NUMLL__ 

C  WRITE  1 21,  115riMTS»NLilH»NPF  ,NUMNP,MJKEL 

C 

__  00  3?-  NTL=!  ,M1S _ _ 

P£Oli*»!2‘'OlNT,TI«u 

C  VRlTEin,  K5C)NT,TTmE 

_ o  o„  *  uuss  1 1  mu  m  N  p _ 

READ C?3, 110" INNODE. f XYZ (I ,N ) ,I=! ,NDTM > , IU 1 1 .N ) , 1  =  1 , NDF I 
continue 

r  VRTTF|?l,i3‘,R) 

WRITE!?;  ,i  1ST)  it  tin,  J)  ,  1  =  1  ,NDF  > , J= 1 , NbPNP ) 

t > 


e:-*? 

r 

r 


Dr*  i,:r 

REA0l?:,iA:.',)»!E'LH,  lYxldll,  J,N)  ,  Jsi.NOlHl  ,1  =  1  ,NGIK) 

CrNT:VUE 

_W  R I  T_^_l  ?  1 , 1 «» J  _ 

THE  H vC°OST A T IC  PRESSURES  ARE  CONVERTED  FROM 
DYNE / I CM**E I  INTO  PST 


«  1*161  l,l,K'l*Xl&t2»2.:,>*XIRlT 
•?03ECu5,Nsi ,NUMEL > 


CONTINUE 


!,.>■»  )/?./ 


ii’C  FO^Htmii) 

U5C..  F09HATUH:!lTS“,Iv»_a2!L».lHNO.lK  =  ,I**,SX,**HNOF  =  ,:U,PX,6HNUKNP:,m,dY, 
.  6H4UKEL=,I4) 

UrZ  FC»PHiTi:E,El?.SL?) 

CI25r  FORKATIcrx.llHTlME  STEP  = , I  5 , 1 ?X , 6MT IKE  =  ,EK.5E?) 

!3EC  F0°H  im?  ,fc£  12.  jp;> 
rl  74R  F0°UATIi.X  .fHPISPL) 

FOPf.ATIor  12. CE?) 

!  4 C u  FORK  A  T 1 1’  ,  6E 1 1  •  SEE  )~ 

CI44R  F0RMATI6X, 17HSTKESS  COMPONENTS) 

1450  F0RKATI6EK.EE?) 

STOP 

END 
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UTL.WRITE80 /ASCII 


l.'TL<!  )  .WP1TE  =  -/Asc:if  :» 

identification  DIVISION. 

PR  06" AM  -  ID  •  ^i%nrASC°.l  . 

. .  "pn'VIPCWMENT'  division."" 

C  CK'FI  Gl'R  ATI  ON  StCTTCM. 

OCUkCL-COMPUTLP.  u N' I V A C -  1 1 " r.  • 

CE JECT-COMPuTtO.  uMI VAC-  1  l"u  .  " 

I N^UT -OUTPUT  srCT;o,\. 
fill-conthol. 

SELECT  InF ILL  ASSIGN  TO  CARD-PEADER. 

SELECT  OUT  ASSIGN  TO  INTERCHANGE  2. 

OATA_DJV!SIwN._. _ _  _ _ _ 

FILE'  SECTION. 

FD  Infill  i.i*2FL  rtFCORD  IS  GKITTFD  RECORDIN' 

CLOCK  CONTAINS _ <«C  RrCOHD 

Cl  "iNF EC  PIC  X(8u). 
r D  Out  l/sel  record  IS  OKITTED 
7£C09CT:.’F  f 
plock  contains  -t?  Records. 

P I  OUTREC  Pir  y<o_}. 

V.'O  "K I  N'3 -STORAGE  SECTION. 

Cl  CONN ’IT  ' ^IC  9**99  ‘  vALUEC. 

PROCEDURE  DIVISION. 

STARTH  .  _ 

CPF  n  INPUT  '  I N FI L E  W IT H '  N 0 ' R El  Nn. 

OPEN  OJ'T °uT  0UT  WITH  NO  PEWIN'L. 

LOOP. 

rE*  D  INFJlr  AT  END  GO  TO  THEVu. 

ADD  *  TO  CGUNTIT. 

'•  O'V  £  IMNEC  to  gutfec. 
i'a:ir  cuTr,rc. . 

G C  TC  LOOP. 

THEN? . 

CLOSE  INFILL  ».  I TH  NC  REWIND. 

CLOSE  Cj  j  V'l  TH  ND  PEW  It'D. 

DISPLAY  *  COUNT  •  C "UNT IT  UPON  PRINTER. 
STOP  Fs,V. 


ID  4/1 


FILKFUP 


14. 

M5 


1:«CI:lETE  E«P1.  »P2.  ,D3. ,D4. ,D5. ,D6.  ,D7.  , D8. , D9. , Die.  ,P11.,D12.  ,DI3.S 

*  .M6.,Pl7.  .MS.  ,P19.  ,P20.{D21. 

2: 8PELETE  E-Pl. ,P2.,P3. ,P4.,P5.,P6, 


14.  r 

p I?. , Fie., pi ?. ,Pie.,pie.,p2e.,P2i. 


4 

5 
€ 

7 

8 
9 

10 

11 

12 

13 

14 

15 
Id 

17 

18 
19 
28 
21 
22 

23 

24 

25 

26 

27 

28 
29 
3G 

31 

32 
73 

?5 

36 

37 
i:i 

3? 

40 

41 

42 

43 

44 

45 

46 

47 

48 
4? 
58 

51 

52 

53 

54 
35 

56 

57 

58 
39 
CO 


•EDIT  F, 
MAK£<1,23)  Dl. 
MAKE(24,8>  PI. 
MAkE<£7,23>  P2. 
MAKE (50,  8)  P2. 
HAKE(53, 23)  D3. 
MAKE<  76, 8)  P3. 
MAKE (79, 23)  D4. 
MAKE (192,8)  P4. 
MAKE <105 
MAKEC128 
MAKE (131 
MAKE( 154 
MAKE< 157 
MAKE <189 
MAKE (1 83 
HAKE <296 
HAKE(209 
MAKE<232 
MAKE <235 
MAKE(2S8 
NaKE (261 
MAKE (284 
MAKE <237 
MAKE (310 
MAKE (313 
MAKE <336 
MAKE (339 
MAKE (362 
MAKE <365 
MAKE (388 
MAKE (391 
MAKE <4 14 
MAKE (41 7 
MAKE *440 
MAKE <443 
MAKE (466 
MAKE <469 
MAKE <492 
MAKE<495 
MAKE (518 
MAKE (321 


23 >  D5. 
8>  P5. 
23)  D6. 
8)  P6. 
23)  D7. 
8)  P7. 
23)  D8. 
8)  P8. 
23)  D9. 
8)  P9. 
23)  DIO. 
8>  Pie. 
23)  Dll. 
8)  Pll. 
23)  D12. 
8)  P12. 
23)  D13. 
8)  P13. 
23)  D14, 
0)  P14. 
£3)  D15. 
8)  F15. 
23)  PIC. 
8)  P16. 
23)  D17. 
8)  PI?. 
23)  D18. 
8>  P18. 
23)  D19. 
8)  P19. 
23)  D20. 
8)  P20. 
23)  D21. 


MAKE(544<  8>  P21. 
EXIT 

•MEHUTILITY 
GEOM  REAP 
ALL. 

N 

FUNC  READ 

PI. 

Y 

WRIT 
Ft. 

READ 
F2. 

WRIT 


, P7. ,P8, ,P9.,P10. ,P11.,P12. ,P13.  • 


61: P2. 
62:  READ 
63: PS. 
64:  v 

65: WRIT 
66:  P3. 
67 S READ 
68:P4. 
69:  Y 
79: WRIT 
71SP4. 
72: READ 
73: P5. 
74J  Y 
75: WRIT 
76:  P5. 
77:  READ 
78: P6. 
79:  Y 

80:mrit 

8i:P6. 
82: READ 
83: P7. 
84!  Y 
85: WRIT 
86: P7. 
67: READ 
88: P8. 
89:  Y 
98: WRIT 
9i:P8. 
92: READ 
93:P9. 
94J  Y 

9S:hrit 

96S P9. 
9?: READ 
98: P10. 

=•? :  y 

ioo: writ 
ioiipio. 

102: REAP 
103:  F  1  1. 
t  a*.:  v 
103: writ 
ioe:pii. 
107: READ 
108:  pi  2. 

!09:v 

lie: writ 
i ii:pi2. 

1 12: READ 
l 13: pit. 
l  U:  v 
1 15: WRIT 
iu:Pi3. 

I  17: READ 
118:P14. 
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119:  y 
120: WRIT 
121:P14. 

122: READ 

123:  Pis. 

124:  Y 
125: WRIT 
12«s:pi5. 

127: READ 
123: PIS. 

129:  Y 
130: WRIT 
13UP16. 

132: READ 

133:P17. 

134:y 

135:«RIT 

136:P17. 

137: READ 
138:P18. 
139:y 
140: writ 
l4l:Pio. 

142: READ 
143: P19. 
144:Y 
145: WRIT 
1461P19. 

147: READ 
148: P20. 
149:Y 
150: WRIT 
13l:P28. 

152: READ 
153:P21. 
154:y 
155: WRIT 
156: P21. 

157: EXIT 
158: 0SECTI0H 
159: ALL. 

168: 

161: ALL. U 
162:21 
163:D1. 
l£4:Dl.U 
165: D2. 
166:02.0 
1 6 1 :  tv  • 

1  6i:  t'3,  u 
169:04. 

170: 1*4.  u 

ir;:t-5. 

172: 05. U 
173:08. 
174:d»-;.u 
1?;-:D7. 

176:  07.  U 
177:DS. 

170:  03.  u 
179:09. 

180:  D9.  U 
181:018. 
182:010. U 
183:D11. 
184:Dll.U 


185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 
19? 
193 

199 

200 
281 
202 

283 

284 

205 

206 
207 
206 
289 
218 
211 
212 

213 

214 

215 

216 

217 

218 
219 
228 
221 
222 

223 

224 

225 

226 

227 

228 
229 
238 

231 

232 

233 

234 
2.35 
236 


240 

241 

242 

243 

244 

245 

246 

247 

248 

249 


D12. 
D12.U 
D13. 
D13.U 
D14. 
D14.U 
D13. 
D15.U 
01 6. 
016. U 
017. 
D17.U 
018. 
018. U 
019. 
012.0 
020. 
D20.U 
D21. 
021.  U 
21 
PI. 
Pl.U 
P2. 
P2.U 
P3. 
P3.U 
P4. 
P4.U 
P5. 
P5.U 
P6. 
P6.U 
P7. 
P7.U 
P8. 
P8.U 
P9. 
P9.U 

Pie. 

pio.u 

Pit. 

Pll.U 

P12. 

P12.U 

P13. 

P13.U 

P14. 

P14.U 

P15. 

P15.U 

PI  6. 

P 16.  U 
P 1 7. 
PI  7.0 
Pi  6. 

Pie.u 

F  19. 

P19.U 

P20. 

F 20.  U 

P21. 

P21.U 
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LOOK 


ffltOVtE 

hI.L.U 

M7.U 

P17.U 

KOTA 

V  -3U  X  25 
COHT 
17  23 
-10  0 
UIEW 
READ 

D13.U 

F18.U 

VIEW 

READ 

D12.U 

P19.U 

VIEW 

READ 

D2G.U 

P2G.U 

VIEW 

READ 

D21.U 

P21.1J 

VIEW 

ANIM 

3 

1  3 

«l 

I 

N 


1 


VIEW 


EXIT 
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1— lea  of  Ttacoelaa ticity 

1*  JPagtic  —  tarlal  -  A  notarial  la  called  claatlc  If  the  current  atreaa 
depend a  on  the  current  atraln  only, 

0(t)  ■  E  *  e(t)  . 

2.  Viscoelastic  aaterlal  -  A  material  la  called  vlacoelaatlc  If  the  current 
atreaa  dependa  on  not  only  the  current  atraln  but  alao  the  entire  past 
deforaation  history, 

«ft 

o(t)  -  /  [e(t  -  t)]  . 

T-0 


3.  The  seat  faelllar  elscoelaatlc  baharlora 
a*  Kelaxatlon 


elastic  solid 

linear  V-E  solid 
tpartlally  relaxed) 


linear  V-B  fluid 
•npletely  relaxed) 

—  ...  -  time 


b.  Creep 
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F  +  T«i-ER(“  +  %3r) 


For  unit  step  displacement  at  t,  u(t)  ■  l(t),  we  can  get  relaxation 
function, 


■  4  *  f1 "  ^)  wt>  ■ 


b.  For  unit  step  loading  at  t,  F(t)  ■  l(t),  we  get  creep  function. 


where 

t0  -  the  time  of  relaxation  of  deflection  under  constant  load, 
Te  -  the  time  of  relaxation  of  load  under  constant  deflection, 
ER  -  relaxed  elastic  modulus,  and 
l(t)  -  unit-step  function,  defined  as 


f  1  when  t  >  0 

l(t)  -  J  1/2  when  t  -  0 

I  0  when  t  <  0 
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To  generalize  the  above  formulation  for  a  simple  bar,  which  actually 
consists  of  Infinite  degrees  of  freedom,  Boltsmann  wrote 


t 

/  c(t 


t)  (*)  dt  • 


t 

/  k(t 


t)  (t)  dT  , 


to  account  for  the  linear  load~def lection  relationship  of  a  simple  bar. 

In  a  complete  description  for  the  linear  viscoelastic  solid,  in  terms  of 
stress  and  strain,  we  can  rewrite 

t  t 

•(t)  -  /  K(t  -  t)  •  da  -  /  K(t  -  r)  •  c  •  dr  ,  (6) 


t  t 

*(t)  *  /  C(t  -  t)  •  do  -  /  C(t  -  t)  •  o  •  dr  . 


5.  Bonllnear  viscoelastic  material  (Green  &  Rlvlln) 


t  t 

«(t)  -  /  ^(t  -  t)  «(t)  dt  +  //  42(t  -  t,  t  -  s)  •  C(T)  •  *(s)  dt  ds 


+  ///  ♦*(*  -  T,  t  -  s,  t  -  r)  •  C(T)  •  *(s)  •  e(r)  dT  ds  dr  . 


e  e  e  a 
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6.  Quaa  1-llner  flicotlMtlc  water! al  (Fung) 


•(t)  ■  /  *(t  -  T)  •  dO 


(e) 


(T) 


-  /  *<t  -  T) 


(#)(#■ 


dt 


(9) 


where 

0  -  stress  at  tine  t, 

■  elastic  response  at  tine  t, 
e  *  strain  weasure, 

•  •  reduced  relaxation  function. 
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APTSNDIX  C.  LIST  OS 


MTA  IBOM  RAH/LOffELACt  ITU 


1.  On  th«  following  experlnenta 

we  have  raw 

data  fron 

chart 

recordera . 

Range 

HOB 

Anlnal 

Date 

Kxploalve 

8K773MHI 

_  (ft) 

W 

Orientation 

092881 

8  lb  TNT 

free  field 

21.3 

3.2 

right  aide  on 

100281 

16  lb  TNT 

free  field 

27.0 

6.0 

right  aide  on 

092981 

32  lb  TNT 

free  field 

44.4 

8.0 

right  aide  on 

093081 

64  lb  TNT 

free  field 

76.3 

10.1 

right  aide  on 

t 

10  fc.  length 

ahock  tube 

58.25 

agalnat  end 

100  grain  prlne  cord 

plate 

120381 

31  ft.  length 

100  grain  prlne  cord 

ahock  tube 

58.25 

120481 

61  ft.  length 

100  grain  prlne  cord 

shock  tube 

58.25 

2.  These  are  experiments  abstracted  fro*  Lovelace  report ,  "Double  Peak  Study" 
on  which  we  have  raw  data. 


Date 

Anlaal 

Number 

Animal 

Loca¬ 

tion 

Explosive 

Range 

(ft) 

Tlae 

Duration 

Between 

Shocks 

(•«> 

Animal 

Orientation 

PART  1  (October  1981) 

102881 

2 

North 

2-8  lb  TNT 

11 

1.7 

right 

aide 

an 

1 

North 

2-8  lb  TNT 

3.6 

right 

aide 

on 

102781 

3 

North 

2-8  lb  TNT 

11 

5.6 

right 

aide 

an 

2 

North 

2-8  lb  TNT 

7.6 

right 

side 

on 

102581 

14 

East 

2-8  lb  TNT 

11 

9.7 

right 

side 

on 

102781 

1 

North 

2-8  lb  TNT 

11 

11.5 

right 

side 

on 

102681 

1 

North 

2-8  lb  TNT 

11 

13.6 

right 

aide 

on 

PART  2  (May  1982) 

051382 

2 

Heat 

8  lb  pentollte  charges 

10.0 

3.8 

right 

side 

on 

1 

East 

8  lb  pentollte  charges 

3.8 

right 

side 

on 

051182 

2 

Nest 

8  lb  pentollte  charges 

10.5 

9.9 

right 

side 

on 

1 

East 

8  lb  pentollte  chargee 

9.8 

right 

side 

on 

051182 

2 

Nest 

8  lb  pentollte  chargee 

10.5 

9.7 

right 

aide 

on 

PART  3  (May  1982) 

051782 

2 

Nest 

8  lb  pentollte  charges 

11 

9.6 

right 

side 

on 

1 

East 

8  lb  pentollte  charges 

9.6 

right 

side 

on 

052182 

2 

Nest 

8  lb  pentollte  charges 

11 

9.6 

right 

side 

on 

1 

East 

8  lb  pentollte  charges 

9.6 

right 

side 

on 

051982 

2 

Nest 

8  lb  pentollte  charges 

11 

3.7 

right 

side 

on 

1 

East 

8  lb  pentollte  charges 

3.7 

right 

side 

on 

052082 

2 

Nest 

8  lb  pentollte  charges 

11 

3.7 

right 

side 

on 

1 

East 

8  lb  pentollte  charges 

3.7 

right 

side 

on 

120 


distribution  list 


12  copies 
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12  copies 


1  copy 


1  copy 
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